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interply  crack  caused  by  interlaminar  stresses.  These  are  the  two  fundamental 
modes  exhibited  in  matrix  cracking. 

Experiments  using  graphite-epoxy  composite  systems  have  established  that 
matrix  cracking  is  actually  a  process  in  which  a  great  number  of  similar  cracks 
occur  during  the  course  of  loading.  The  multiplicity  in  the  crack  formation 
process  stems  from  the  existence  of  material  flaws  which  distribute  randomly 
throughout  the  laminate,  and  from  the  internal  reinforcement  structure  which 
inhibits  individual  crack  from  growing  into  large  proportions. 

When  viewed  at  the  ply  level,  growth  of  a  matrix  crack  is  essentially 
brittle  in  nature;  its  growth  path  initially  follows  one  of  the  two  fundamental 
modes.  When  blunted,  the  initial  mode  may  turn  to  the  other  mode.  Thus,  the 
individual  growth  path  may  be  considered  as  piece-wise  self-similar. 

The  general  methodology  developed  in  this  study  invokes  the  fundamental 
concepts  of  ply  elasticity  which  assumes  that  the  material  ply  in  the  laminate 
is  individually  elastic,  homogeneous  and  anisotropic.  The  constitutive  proper¬ 
ties  can  be  represented  by  a  set  of  effective  constants  and  the  failure  proper¬ 
ties  are  governed  by  randomly  distributed  material  flaws.  The  latter  is  char¬ 
acterized  by  an  effective  flaw  distribution  which  is  considered  a  basic  ply 
property.  Thus,  for  each  individual  flaw,  its  propagation  behavior  is  governed 
by  criteria  derived  from  the  classical  fracture  mechanics,  while  the  flaw 
distribution  enables  the  description  of  the  multiplicity  and  the  stochastic 
nature  of  the  matrix  cracking  process. 

This  report  highlights  these  methodical  developments,  along  with  illustra¬ 
tive  experimental  and  analytical  results. 

Three  technical  papers  are  appended  to  this  report.  These  papers  contain 
the  details  of  the  development  and  much  more  results.  Although  each  paper  was 
written  in  self-contained  format,  they  collectively  represent  the  bulk  of  the 
technical  development  accomplished  in  this  research  program. 
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ABSTRACT 


^  This  report  highlights  the  development  o£  a  unified  method  for  the  formation 
mechanisms  of  matrix-dominated  cracks  in  a  class  of  structural  laminates  that  are 
made  using  unidirectionally  fiber  reinforced  composite  plies. 

Because  of  the  collimated  reinforcement  structure  exhibited  in  this  class  of 
laminates,  the  growth  of  an  individual  matrix  crack  is  confined  either  in  the 
fiber-matrix  interface  or  in  the  ply-to-ply  interface.  The  former  occurs  as 
intraply  crack  caused  primarily  by  in-ply  stresses,  while  the  latter  occurs  as 
interply  crack  caused  by  interlaminar  stresses.  These  are  the  two  fundamental 
modes  exhibited  in  matrix  cracking. 

Experiments  using  graphite-epoxy  composite  systems  have  established  that 
matrix  cracking  is  actually  a  process  in  which  a  great  number  of  similar  cracks 
occur  during  the  course  of  loading.  The  multiplicity  in  the  crack  formation 
process  stems  from  the  existence  of  material  flaws  which  distribute  randomly 
throughout  the  laminate,  and  from  the  internal  reinforcement  structure  which 


inhibits  individual  crack  from  growing  into  large  proportions. 


’  When  viewed  at  the  ply  level,  growth  of  a  matrix  crack  Is  essentially  brittle 
in  nature;  its  growth  path  initially  follows  one  of  the  two  fundamental  modes. 

When  blunted,  the  Initial  mode  may  turn  to  the  other  mode.  Thus,  the  individual 
growth  path  may  be  considered  as  piece-wise  self-similar.^ _ ^ 

The  general  methodology  developed  in  this  study  invokes  the  fundamental 
concepts  of  ply  elasticity  which  assumes  that  the  material  ply  in  the  laminate 
is  Individually  elastic,  homogeneous  and  anisotropic.  The  constitutive  properties 
can  be  represented  by  a  set  of  effective  constants  and  the  failure  properties  are 
governed  by  randomly  distributed  material  flaws.  The  latter  is  characterized  by 
an  effective  flaw  distribution  which  is  considered  a  basic  ply  property.  Thus, 
for  each  individual  flaw,  its  propagation  behavior  is  governed  by  criteria 
derived  from  the  classical  fracture  mechanics,  while  the  flaw  distribution  enables 
the  description  of  the  multiplicity  and  the  stochastic  nature  of  the  matrix 
cracking  process. 

This  report  highlights  these  methodical  developments,  along  with  illustrative 
experimental  and  analytical  results. 

Three  technical  papers  are  appended  to  this  report.  These  papers  contain 
the  details  of  the  development  and  much  more  results.  Although  each  paper  was 
written  in  self-contained  format,  they  collectively  represent  the  bulk  of  the 
technical  development  accomplished  in  this  research  program. 


Foreword 


This  is  the  final  technical  report  for  Phaae-II  Study  under  the 
Contract  No.  F49620-79-C-0206  froa  the  Air  Force  Office  of  Scientific 
Research.  The  contract  was  led  to  Drexel  University  with  a  sub-contract 
to  Lockheed  Palo  Alto  Research  Laboratory,  covering  the  period  from 
September  1,  1979  to  September  30,  1983. 

The  contracted  research  program  consisted  two  phases.  The  Phase-I 
Study  spanned  the  period  from  September  1,  1979  to  September  30,  1981;  and 
the  Phase-II  Study  spanned  the  period  from  October  1,  1981  to  September  30 
1983. 

Or.  A.  S.  0.  Vang  of  Drexel  University  was  the  principal  investigator 
and  Dr.  F.  V.  Crossman  of  Lockheed  was  the  co-principal  investigator. 

Dr.  Wang  was  assisted  by  his  former  and  present  students,  principally 
Dr.  G.  E.  Law,  Jr.,  Dr.  N.  N.  Kishore  and  Mr.  C.  A.  Li.  Dr.  Crossman  was 
assisted  by  Mr.  V.  J.  Warren. 

The  authors  would  like  to  thank  Dr.  Anthony  K.  Amos  and  Major  David  A 
Glasgow,  both  of  AFOSR.  Dr.  Amos  waa  the  contract  monitor  for  Phaae-I 
Study  and  Major  Glasgow  was  the  contract  monitor  for  Phase-II  Study. 


I.  INTRODUCTION 


This  report  contains  the  results  of  a  research  program  conducted  under 
the  auspices  of  the  Air  Force  Office  of  Scientific  Research.  The  stated 

< 

objectives  of  the  research  are  (1)  to  investigate  the  basic  mechanisms  of  sub¬ 
laminate  cracks  in  a  class  of  aerospace  structural  laminates;  and  (2)  to 
develop  a  rational  modelling  procedure  by  which  the  formation  and  growth 
processes  of  the  observed  sublaminate  cracks  can  be  described  in  terms  of 
engineering  quantities. 

The  particular  class  of  laminates  considered  in  this  study  included  those 
fabricated  using  continuous  fiber  reinforced  resin  Bystems,  e.g.  the 
unidirectional  graphite-epoxy  prepreg  tape.  Sublaminate  failure  in  this  class 
of  laminates  begins  generally  with  local  matrix  cracks  at  low  stress  levels. 
The  individual  growth  paths  of  these  cracks  almost  always  follow  either  the 
fiber/matrix  interface  or  the  ply-to-ply  interface.  Hence,  there  are  only  two 
fundamental  modes  for  matrix  crack  growth,  the  intraply  cracking  mode  and  the 
interply  cracking  mode.  The  present  research  is  concerned  with  the  physical 
mechanisms  of  the  two  fundamental  cracking  modes  as  well  as  their  mutual 
interactions. 

Consequently,  a  coup  rehens ive  effort  was  made  during  the  research  period 
which  embodied  both  an  experimental  and  an  analytical  investigation. 

The  first  phase  of  the  research  was  directed  to  the  basic  mechanics  of 
intraply  and  interply  crackings  individually,  whereby  a  unified  energy  method 
was  developed  for  their  formation  mechanisms.  Several  technical  papers  were 
published  on  one  or  more  aspects  related  to  the  development  of  the  analytical 
model  (see  Refs.  (1]  to  [8],  Appendix  A).  An  overall  exposition  of  the  first 
phase  study  results  was  contained  in  Ref.  {9],  Appendix  A. 


The  second  phase  of  the  research  was  directed  to  the  understanding  of 
multiple  cracks  and  cracking  modes  interaction*  These  sublaminate  failures 
occur  following  the  onset  of  one  of  the  fundamental  modes  of  matrix  cracking. 

Multiple  cracks  and  their  mutual  interactions  presented  both  conceptual  ■ 
and  physical  difficulties.  The  former  is  related  to  the  rational  derivation 
of  the  desired  analytical  model,  while  the  latter  is  related  to  the 
mathematical  and/or  numerical  techniques  necessary  to  execute  the  model  for 
real  problems.  In  this  regard,  several  major  developments  are  noted.  First, 
the  concept  of  "effective  flaw  distribution"  was  introduced  which  replaces  the 
conventional  strength  concept  for  the  basic  ply  material.  It  is  assumed  that 
these  flaws  cause  multiple  sublaminate  cracks.  Rational  arguments  based  on 
experimental  observations  and  simple  testing  techniques  are  devised  in  order 
to  characterize  these  flaws  as  a  ply  material  property. 

Secondly,  computer  simulation  procedures  based  on  3-D  finite  element 
method  and  incorporating  fracture  mechanics  criteria  were  developed  in  order 
to  describe  the  interacting  and  stochastic  nature  of  nultiple  cracks. 

Finally,  extensive  experiments  were'  conducted  in  concert  with  each  of  the 
analytical  development  phases.  These  experiments  involve  the  most  up-to-date 
destructive  and  non-destructive  damage  inspection  techniques,  which  provide 
not  only  a  visual  description  of  the  various  matrix  cracks  in  the  laminates 


but  also  a  chronological  documentation  of  the  crack  formation  processes  during 
the  course  of  loading.  Both  are  essential  to  the  initial  formulation  and  the 
final  correlation  of  the  analytical  model. 


This  report  contains  the  results  of  the  second  phase  of  research. 


A’  N*  %\ 


II.  RESULTS 


Major  results  obtained  in  the  second  phase  study  have  been  presented  in  a 
series  of  three  technical  papers.  These  papers  are  listed  in  Appendix  A  as 

I 

Ref.  [10],  [11]  and  [12],  and  are  appended  to  this  report  as  Appendix  B,  C  and 
D. 

The  first  paper  (in  Appendix  B)  is  entitled  "Fracture  Mechanics  of  Sub- 
Laminate  Cracks  in  Composite  Laminates."  In  this  paper,  a  detailed  account  is 
given  for  the  mechanics  of  transverse  cracking  and  for  the  mechanics  of  free 
edge  delamination;  the  concept  of  effective  flaws  is  augmented  to  simulate 
multiple  cracks  and  their  mutual  interactions.  Results  from  two  case  studies 
are  included  to  provide  a  comparison  between  the  analysis  and  the  experiment. 
In  the  first  case  study,  90° -ply  transverse  cracks  in  the  family  of  [02/90n]s, 
n  *  1,  2,  3,  laminates  are  tested  and  sinulated.  In  the  second  case  study, 
the  family  of  [+25/90n]s,  n  *  1/2,  1,  2,  3,  4,  6,  8,  laminates  are  tested  and 
simulated.  It  is  noted  that  de lamination  mode  is  not  involved  in  the  first 
case  study,  while  both  transverse  cracking  and  delamdLnatlon  are  Involved  in 
the  second  case  study.  Depending  on  the  value  of  n  of  the  latter  case,  one 
cracking  mode  can  be  induced  before  the  other.  Generally,  the  analytical 
model  can  reliably  predict  the  various  cracking  events. 

The  second  paper  (in  Appendix  C)  is  entitled  "Influence  of  Ply  Thickness 
on  Damage  Accumulation  and  Final  Fracture. ”  The  paper  documents  the 
experimental  finding  for  a  series  of  graphite-epoxy  laminate  families  in  the 
form  l^/^Onlg,  where  the  thickness  of  the  90°-ply  varies  from  2  to  16  plies 
(a  *  1  to  8)  and  the  angle  "6"  takes  on  values  of  0®,  25°,  30°  and  60®. 
Incremental  tensile  loading  followed  by  x-radiographic  damage  inspection  were 
conducted  to  document  chronologically  damage  accumulation  in  the  laminates; 
deplying  experiments  and  real-time  x-ray  during  tensile  loading  were  also  used 


to  monitor  the  details  nearing  imminent  laminate  failure.  It  is  shown  in  this 
experimental  study  that  matrix  cracking  modes  and  their  propagative  behavior 
change  with  laminate  stacking  sequence  and  ply  thickness  in  a  manner  that 
supports  the  analytical  model  based  on  the  criteria  of  fracture  mechanics.  1 
Furthermore,  the  effects  of  sub-laminate  cracks  on  the  primary  load-carrying 
plies  were  also  evaluated  experimentally. 

The  final  paper  (in  Appendix  D)  is  entitled  "On  Crack  Development  in 
Graphite  Epoxy  [02/90nl8  Laminates  Under  Uniaxial  Tension-  "  This  paper  makes 
use  of  the  experimental  data  reported  in  the  previous  paper.  The  particular 
laminates  chosen  belong  to  the  family  of  [02/90n]s,  n  =  1,  2,  and  4.  When 
these  laminates  are  subjected  to  uniaxial  tension,  transverse  cracks  in  the 
90°-ply  develop  at  tensile  strain  as  low  as  0.  3Z.  As  the  applied  tension 
increases,  the  cracks  continue  to  form  in  increasing  numbers.  At  some 
critical  late  stage  of  loading,  however,  several  other  modes  of  matrix 
cracking  emerge,  all  stemming  from  the  sites  of  the  transverse  cracks.  Among 
them  are  (1)  splitting  cracks  in  the  0°-ply  that  span  across  one  or  more 
90°-ply  transverse  cracks,  (2)  localized  delamination  in  the  0°/90°  interface, 
emanating  from  the  point  where  a  0°-ply  split  intersects  a  90°-ply  transverse 
crack;  and  (3)  localized  edge  delamination  in  the  0e/90°  interface  emanating 
from  the  corner  point  where  a  90° -ply  transverse  crack  meets  the  free  edge  of 
the  laminate  coupon.  These  late  stage  matrix  crackings  seem  to  precipitate 
weakening  of  the  load-carrying  ply  of  the  laminate,  resulting  in  lowered 
apf  'ent  strength.  This  paper  attempts  to  analytically  describe  all  these 
Important  cracking  events  by  means  of  the  previously  developed  fracture 
mechanics/ef fective  flaw  method.  Emphasis  is  placed  on  explaining 
qualitatively  and  quantitatively  why  a  certain  mode  of  cracking  occurs  while 
other  modes  do  not;  especially  for  the  cracking  state  which  exists  during  late 


stage  loading,  an  analysis  is  made  to  explain  why  that  state  is  critical  to 
the  apparent  laminate  strength  reduction.  As  stated,  all  predictions  are 
compared  to  the  experimentally  obtained  data. 

Although  each  of  the  three  separate  papers  is  written  in  self-contained  • 
format,  they  do  maintain  the  essential  continuity  in  contents  and  consistency 
in  presentation.  Collectively,  these  three  papers  represent  the  bulk  of  the 
second  phase  study  results. 

A  list  of  technical  presentations  made  during  both  the  first  phase  and 
the  second  phase  studies  is  included  in  Appendix  A. 
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ABSTRACT 

This  paper  presents  an  overview  of  a  fracture  mechanics  approach  to 
some  of  the  most  frequently  encountered  matrix-dominated,  sub-laminate 
cracks  in  epoxy-based  composite  laminates.  By  "sub-laminate",  it  is  meant 
that  the  cracks  are  internal  to  the  laminate,  generally  invisible  macro- 
scopically;  but  are  much  larger  in  size  than  those  microcracks  considered 
in  the  realm  of  micromechanics.  The  origin  of  sub-laminate  cracks  is 
assumed  to  stem  from  the  coalesence  of  natural  material  flaws  (also  micro¬ 
cracks)  which  occur  under  a  certain  favorable  laminate  stress  condition. 

a 

Thus,  the  modeling  of  the  mechanisms  of  sub-laminate  crack  Initiation  and 
propagation  is  essentially  mechanistic  and  probabilistic  in  nature.  Some 
specific  results  from  several  analytical/experimental  investigations  using 
graphite-epoxy  laminates  are  presented  and  discussed  in  this  paper. 

KEY  WORDS,  macroscopic  analysis;  fracture  mechanics;  transverse  cracks; 
delamination;  mechanistic  modeling;  stochastic  simulations;  sub-laminate 
cracks  interaction;  final  failure  mechanisms;  failure  modes  sequence. 


INTRODUCTION 


Failure  mechanisms  In  structural  composite  laminates  have  been  viewed 
at  several  dimensional  levels.  Consider,  for  instance,  the  curved  lamina¬ 
ted  panel  shown  In  Figure  1(a).  Failure  of  this  structural  component  may 
be  caused  by  a  loss  of  the  global  stiffness  when  the  applied  load  reaches 
a  certain  critical  value.  To  describe  analytically  the  associated  failure 
mechanisms  and  hence  to  determine  the  critical  load,  structural  mechanics 
methods  such  as  buckling  and  post-buckling  theories  are  employed,  which 
relate  the  change  of  the  structural  geometry  to  the  applied  load.  In  this 
type  of  analysis,  the  stiffness  of  the  laminated  panel  and  the  kinematics 
of  its  displacements  are  among  the  most  predominant  factors. 

On  the  other  hand,  the  same  panel  may  fail  due  to  a  loss  of  strength 
at  a  local  defect;  for  example,  at  a  bolt-hole.  Figure  1(b).  In  this  case, 
rupture  of  material  will  begin  at  the  hole  and  may  propagate  into  a  large 
crack  whenever  a  certain  load  condition  is  reached.  But  in  order  to  define 
the  conditions  for  rupture  initiation  and  propagation,  a  knowledge  of  the 
actual  stress  field  around  the  hole,  and  the  physical  mechanisms  of  the 
material  rupture  process  is  essential.  To  this  end,  it  will  be  necessary 
to  focus  further  on  the  local  defect  area  at  a  much  smaller  dimensional 
level.  As  illustrated  in  Figure  1(c),  for  instance,  the  lamination  de¬ 
tails  of  the  panel,  such  as  layer  interfaces  and  the  stacking  sequence 
are  now  identifiable.  Consequently,  quantities  of  size  comparable  to  the 
layer  thickness  become  important;  the  influence  of  a  local  interface  de¬ 
fect,  a  transverse  crack,  a  fiber  split  in  the  0-layer,  etc.  are  new  fac¬ 
tors  to  be  considered.  For  it  is  believed  that  the  development  of  rupture 
around  the  hole  is  precipitated  by  these  sub-laminate  cracking  activities. 


Any  failure  analysis  performed  at  this  dimensional  level  must  address 
the  individual  mechanisms  of  the  various  sub-laminate  cracks  as  well  tea 
their  Interaction  and  coalescence  mechanisms  as  they  occur  simultaneously 
and/or  sequentially. 

It  is  well  known  that  the  formation  and  propagation  mechanisms  of  a 
crack  are  governed  by  the  conditions  that  exist  in  the  close  vicinity  of 
the  crack  front.  In  the  case  of  the  sub-laminate  cracks  mentioned  above, 
further  focusing  of  the  crack  front  region  will  reveal  the  microstructure 
of  the  composite  system.  Figure  1(d).  Here,  a  clear  distinction  can  be 
made  between  the  fiber  and  the  matrix  phases  of  the  material.  At  this 
dimensional  scale,  one  also  finds  a  random  distribution  of  material  micro¬ 
flaws,  be  it  voids,  broken  fibers,  matrix  crack,  debonded  fiber-matrix  in¬ 
terfaces,  just  to  mention  a  few,  see  Figure  1(e). 

Although  the  Inherent  mlcroflaws  are  small  in  size,  usually  on  the 
order  of  the  fiber  diameter,  they  may  behave  as  local  stress  risers  under 
the  applied  load.  When  a  certain  local  condition  is  reached,  the  micro¬ 
flaws  interact  and  coalesce  to  form  an  actual  crack,  much  larger  in  size 
and  identifiable  at  the  sub-laminate  level.  Hence,  a  proper  analysis  of 
the  mechanisms  must  be  capable  of  delineating  the  individual  behavior  of 
the  various  types  of  stress  risers,  and  their  Interactions  when  one  is  lo¬ 
cated  near  another.  All  these  will  depend  profoundly  on  the  probabilistic 
nature  of  the  mlcroflaws  in  terms  of  their  size,  location  and  density 
distributions. 

The  physical  process  of  material  failure  as  portrayed  in  the  preced¬ 
ing  discussions  is  seen  to  span  a  wide  range  of  dimensions.  Though  the 
entire  process  is  essentially  one  continuous  event,  failure  analysis  could 


be  conducted  only  within  isolated  dimensional  regimes.  As  illustrated  in 
Figure  1,  the  entire  dimensional  spectrum  is  separated  roughly  into  three 
analysis  regimes;  namely  the  structural  mechanics,  the  macromechanics  and 
the  micromechanics  regimes.  Within  each  analysis  regime,  the  material 
failure  process  involves  a  distinct  set  of  material  and  geometrical  para¬ 
meters. 

Analytical  and  experimental  treatments  on  laminate  failure  at  the 
structural  mechanics  level  have  generally  been  approached  by  using  the 
classical  lamination  theory  (see  (1]).  Stiffness-predominant  structural 
responses  such  as  panel  buckling  and  post-buckling  have  been  successfully 
analyzed  [2].  The  effect  of  local  material  damage  (e.g.  from  a  bolt-hole, 
a  delamination  caused  by  impact,  etc.)  upon  the  global  response  can  also 
be  evaluated  [3].  But,  a  more  detailed  description  about  the  damage  form¬ 
ation  and  propagation  mechanisms  always  required  separate  treatment. 

For  a  class  of  epoxy-baaed  composite  laminates,  e.g.  graphite-epoxy 
systems,  material  damage  mechanisms  around  a  bolt-hole  or  near  the  tip  of 
a  line-notch  have  been  modeled  by  an  empirical  formulae  which  is  de¬ 
rived  from  the  classical  fracture  mechanics  for  structural  metals  [4-6]. 
Essentially,  the  notched  laminate  is  regarded  as  a  2-dimensional  body, 
and  failure  from  the  notch  is  assumed  in  the  form  of  a  line  crack  propaga¬ 
tion.  In  order  for  the  fracture  formula  to  correlate  with  a  body  of  ex¬ 
perimental  data,  the  notch  size  parameter  must  be  adjusted  to  include  an 
empirically  defined  "intense  energy  region"  at  the  crack  tip.  Thus,  to¬ 
gether  with  the  laminate's  fracture  toughness  (Kc  or  Gc) ,  yet  to  be  fit¬ 
ted  experimentally,  the  model  consists  of  two  disposable  parameters 
which  can.  be  adjusted  only  when  a  suitable  set.  of  test  data  is  given. 


Clearly(  the  approach  Is  totally  empirical  in  nature.  The  tern 
"intense  energy  region"  ie  in  fact  a  gross  representation  of  local  sub¬ 
laminate  cracking  activities  at  the  ply  level.  Consequently,  both  the 
"intense  energy  region"  and  the  "fracture  toughness"  parameters  may  differ 
from  one  lamination  to  another  lamination.  It  would  be  ideal,  therefore, 
to  be  able  to  perform  failure  analysis  at  the  ply  level  so  as  to  understand 
how  sub-laminate  cracks  occur,  and  how  they  lead  to  damage  development  near 
the  bolt-hole  or  the  line-notch  in  the  laminate. 

Early  studies  on  failure  in  fiber  composite  systems  have  dealt  mainly 
mechanisms  in  the  micromechanics  regime.  For  instance,  a  considerable  amount 
of  theoretical  and  experimental  treatments  was  given  to  the  subject  of  fiber 
matrix  Interface  mechanics  C7  -lo].  Fracture  models  for  such  micro flaws  as 
those  depicted  in  Figure  1(e)  have  been  characterized  in  terms  of  energy 
absorption  processes  Cll3.  Although  modeling  the  various  types  of  microflaws 
is  essential  for  an  understanding  of  the  strengthening  mechanics  of  the 
composite  system,  these  "micro"  models  have  not  yet  provided  a  further 
analysis  on  the  mutual  interaction  and  the  coalesoenoe  mechanisms  amongst 
the  microflaws.  Indeed,  it  would  be  ideal  again  to  have  a  unified  model 
which  not  only  accounts  for  the  mechanistic  failure  behavior  but  also  the 
statistic  nature  of  the  microflaws.  Such  an  analysis  can  then  bring  the 
knowledge  about  microcrack  mechanisms  to  the  regime  of  macromechanics. 

Several  mechanistic/probabilistic  failure  models  have  been  developed 
in  the  past.  Zweben  and  Rosen  £123  studied,  for  example,  the  tensile 
strength  properties  of  unidirectional  composites  by  assuming  the  micro¬ 
flaws  in  the  form  of  local  broken  fibers.  Though  the  basic  concept  of 
material  flaw  distribution  can  theoretically  be  extended  to  multi- 


directional  laminated  systems,  the  complexity  in  geometry  as  well  as  the 
multiplicity  In  the  sub-laminate  cracking  modes  make  this  type  of  approach 
so  far  unproductive. 

In  recent  years,  increased  attention  has  been  given  to  fracture  model¬ 
ing  within  the  frame  work  of  ply-elasticity  and  the  theory  of  the  classical 
fracture  mechanics.  At  this  level  of  analysis,  the  individual  material  plies 
in  the  laminate  are  replaced  by  homogeneous  and  anisotropic  bodies,  where 
no  distinction  of  the  fiber  and  the  matrix  phases  need  be  made.  Specifically, 
the  stiffness  property  of  the  material  ply  is  represented  by  a  set  of  "effec¬ 
tive"  constants  to  be  determined  experimentally |  while  the  strength  property 
is  governed  by  a  distribution  of  material  flaws  whose  fracture  events  lead 
to  multiple  failures  at  the  sub-laminate  level.  Since  by  a  "homogeneous" 
approximation  the  true  identities  of  the  inherent  material  flaws  are  lost, 
their  effects  on  the  formation  and  propagation  of  multiple  cracks  are  re¬ 
tained  through  an  introduction  of  a  new  concept*  namely,  the  concept  of 
"effective"  material  flaw  distribution,  as  a  basic  ply  property. 

As  it  will  be  discussed  in  more  detail  later  in  this  paper,  the  macro¬ 
scopic  view  point  incorporating  the  effective  flaw  concept  enables  a  rational 
formulation  of  a  unified  fracture  model  which  is  applicable  to  a  class  of 
sub- laminate  cracking  phenomena  observed  in  brittle  matrix  systems,  such  as 
graphite-epoxy  composites. 

Early  investigations  on  composite  failure,  notably  by  Corten  [13],  Wu 
[14],  Kanninen,  et.  al.  [15],  have  already  articulated  the  viability  of  such 
an  approach.  Advances  have  since  been  made,  due  largely  to  the  ever-expanding 
computational  capabilities  and  the  ever-revealing  damage  inspection  tech¬ 
niques.  These  modem  facilities  have  provided  a  means  for  more  rational 
correlations  between  experiments  and  analyses. 


In  a  series  of  papers  [16  -  22],  Wang,  Crossman,  et.  al.  developed  a 
unified  energy  model  within  the  context  of  macromechanics.  It  describes  the 


growth  mechanisms  in  a  class  of  matrix-predominant  sub-laminate  cracks. 

The  specific  cracks  considered  included  interlaminar  (delamination)  and 
intralaminar( transverse)  crackings.  These  sub-laminate  failure  modes  have 
been  found  most  prevalent  in  graphite-epoxy  laminates. 

The  main  objective  of  this  paper  is  to  give  an  overview  on  the  theore¬ 
tical  development,  and  to  present  some  of  their  major  results. 

FRACTURE  CONSIDERATIONS  IN  MACROMECHANICS 
The  energy  model  developed  by  Wang  and  Crossman  is  formulated  on  the 
energy  release  rate  concept  of  the  linear  elastic  fracture  mechanics 
(LEFM) .  All  variables  in  the  model  are  defined  within  the  context  of  the 
macromechanics.  Since  fibrous  composites  possess  some  unique  characteris¬ 
tics,  special  considerations  had  to  be  given  to  defining  some  of  the  varia¬ 
bles. 

Application  of  the  Griffith  Criterion 

The  essence  of  the  LEFM  is  that  material  failure  is  not  defined  from 
the  stand-point  of  strength  as  a  constant  material  property.  Rather,  it 
is  determined  by  an  analysis  of  the  kinetics  of  the  actual  process  of 
fracture  propagation.  The  classical  result  of  Griffith  [23]  pertained  to 
an  elastic  plate  which  is  uniformly  stretched  in  one  direction  by  o.  The 
plate  has  a  through  crack  of  size  2a,  orientated  normal  to  the  direction 
of  o.  The  length  of  the  crack  is  assumed  small  but  finite;  and,  the 
material  is  homogeneous  and  isotropic.  Griffith  postulated  that  at  the 
instance  of  crack  extension,  a  loss  of  the  stored  elastic  strain  energy  U 
results;  and  this  energy  is  assumed  to  be  converted  entirely  into  the 


surface  energy  of  the  crack.  A  balance  of  energy  during  a  stable  crack 
extension  leads  to  the  criterion 


1?.  .  ii 
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Equation  (1)  defines  the  general  condition  under  which  the  existing 
crack  begins  self-similar  propagation  in  the  direction  of  the  crack. 

For  the  plate  problem,  the  critical  stress  at  the  instance  of  crack 
extension  can  be  derived  from  (1) ,  giving 


a 
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where  y  is  the  free  surface  energy  density  of  the  material. 

When  applied  to  engineering  problems,  the  Griffith  theory  is  often 
modified  for  practical  considerations.  For  instance.  Equation  (2)  becomes 
unbounded  as  a  0.  Certainly,  no  real  material  can  sustain  an  infinite 

stress.  This  limitation,  however,  can  be  circumvented  by  invoking  the 

existence  of  material  flaws.  That  is  to  assume  for  the  material  some 

characteristic  distribution  of  flaws;  the  worst  of  which,  having  a  size 

of  2a  *  acts  like  a  real  crack.  It  then  determines  a  finite  critical 
o 

stress  according  to  Equation  (2). 

Of  course,  flaws  do  exist  in  real  materials,  especially  in  fibrous 
composite  systems.  But  the  physical  identity  is  lost  at  the  dimensional 
level  where  the  analysis  is  performed.  The  quantity  aQ  can  be  defined 
only  empirically  as  an  effective  material  property.  As  it  will  be  dis¬ 
cussed  later,  the  Quantity  &Q  Is  in  fact  a  random  variable,  representing 
the  size  distribution  of  the  material  flaws. 
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Another  practical  consideration  is  related  to  the  definition  of  y, 
the  free  surface  energy  density  of  the  material.  For  crack  in  brittle 
material  such  as  glass,  then  y  is  as  defined.  For  most  other  engineering 
materials,  crack  extension  is  found  to  associate  a  certain  degree  of  in¬ 
elastic  deformation  near  the  crack-tip  region.  Furthermore,  the  crack 
extension  path,  or  the  crack  surface,  shows  a  certain  degree  of  rugged¬ 
ness,  depending  on  the  heterogeneity  of  the  material  viewed  at  a  micro¬ 
scale.  Early  studies  by  Irwin  [24],  Orowan  [25]  and  others  on  structural 
metals  considered  the  right-hand  side  of  Equation  (1)  the  irreversible 
work  required  to  create  a  unit  crack  surface  area.  Hence,  the  quantity  y 
can  be  interpreted  as  the  energy  dissipated  in  the  crack-tip  region  during 
crack  extension.  Clearly,  y  will  then  depend  on  the  Inelasticity  as  well 
as  the  microscopic  heterogeneity  of  the  material  locally  (near  the  crack- 
tip).  It,  therefore,  must  also  be  regarded  as  an  Intrinsic  property  of 
the  material  to  be  defined  at  the  macromechanics  level.  Conceivably,  y 
has  to  be  an  averaged  value  over  a  relatively  large  crack  surface  area 
for  fibrous  composites  in  order  for  it  to  be  a  macroscopic  property. 

These  considerations  are  of  fundamental  importance  when  a  crack-like 
failure  is  modeled  at  the  macroscopic  dimensional  level.  For  only  in  this 
context  can  a  macroscopic  fracture  model  be  developed  along  the  rational 
arguments  of  the  classical  fracture  mechanics. 

In  common  practice,  the  right-hand  side  of  Equation  (1)  is  replaced 
by  Gc  (-  2y),  known  as  the  critical  energy  release  rate  of  the  material. 
The  left-hand  side  is  a  function  of  the  applied  load,  the  geometry  of  the 
body  and  the  size  of  the  crack.  Thus,  for  a  crack  propagating  in  its  own 
direction,  the  Griffith  criterion  is  expressed  commonly  as, 


(3) 


G(<*,a)  =  Gc 

In  view  of  the  foregoing  discussions,  a  successful  model  for  some  of 
the  sub-laminate  cracks  may  be  formulated  if  a  judicious  use  of  Equation  (3) 
is  made.  But,  unlike  in  the  classical  fracture  problems  where  the  crack  length 
a  and  the  fracture  toughness  Gc  are  usually  well  defined,  these  two  Quantities 
are  generally  difficult  to  determine  for  fibrous  composites.  In  addition, 
calculation  of  the  function  G(d,a)  for  a  crack  propagating  in  layered  media 
can  also  be  mathematically  complicated}  invariably,  some  numerical  procedure 
is  required.  All  these  aspects  will  be  further  elaborated  in  the  following. 

The  Concept  of  Effective  Flaw  Distribution. 

In  laminate  stress  analysis,  the  fundamental  approach  has  been  the  so- 
called  ply-elasticity  theory  £l].  The  theory  pertains  primarily  to  laminates 
that  are  composed  of  unidirectional  plies  pressed  together  with  each  ply 
orientated  in  some  designed  direction.  The  individual  plies,  though  con¬ 
sisting  fiber  and  matrix  phases,  are  assumed  macroscoplcally  homogeneous 
having  uniform  properties.  This  assumption  enables  characterization  of  the 
unidirectional  ply  by  two  sets  of  material  constants.  The  first  set  consists 
of  the  ply  deformation  property}  and  the  second  set  is  associated  with  the 
ply  strength  property.  The  stiffness  matrix  in  the  generalized  Hooke's 
law,  and  the  strength  constants  in  the  Quadratic  failure  criterion  are  prac¬ 
tical  examples t 


tf.  =  C.  .  e. 

i  ij  J 


i,J  =  1,6 
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It  is  noted  that  the  constants  C^,  F^,  etc.  are  "uniform'  properties 
of  the  ply  material.  And,  In  particular,  the  Quadratic  failure*  criterion 
governs  only  the  first  failure  of  the  ply.  But,  in  a  laminate,  multiple 
failures  in  a  single  ply  can  occur  at  a  range  of  ply  stress  magnitudes. 

Thus,  the  ply  failure  property  cannot  be  assumed  uniform. 

The  conceptual  argument  for  sub-laminate  cracks  is  that  randomly  dis¬ 
tributed  material  flaws  exist  in  the  microstructures  of  the  unidirectional 
ply;  the  flaw  distributional  characteristics  are  Inherent  to  the  basic 
material  system  used  and,  perhaps  also  to  the  laminate’s  fabrication  pro¬ 
cesses.  But,  within  the  frame  work  of  ply-elasticity,  the  true  Identity  of 
the  inherent  material  flaws  cannot  be  practically  determined.  The  alternative 
is,  therefore,  the  concept  of  "effective"  flaws. 

This  concept  assumes  that  a  characteristic  distribution  of  "effective" 
flaws  can  be  determined  for  a  given  representative  volume  of  the  unidirectional 
ply;  each  of  the  flaws  can  act  like  a  small  crack  in  the  same  sense  as  the 
Griffith  crack.  So,  individually,  the  behavior  of  the  flaw  can  be  described 
by  the  criterion  (3).  Collectively,  they  form  the  source  of  the  observed 
sub-laminate  cracks. 

To  define  the  effective  flaw  distribution,  let  a  be  the  size  of  the  flaw, 
and  S  be  the  relative  position  between  two  neighboring  flaws.  Both  a  and  S 
are  generally  random  variables.  The  probability. density .functions  f( a)  and 
f(s)  then  define  completely  the  effective  flaw  distribution. 

For  each  effective  flaw,  the  size  a,  the  geometry  and  location  relative 
to  the  layering  structure  of  the  laminate  are  defined.  It  remains  then  to 
suitably  calculate  the  available  energy  release  rate  function  G(tf,a)  under 
the  given  laminate  load  d.  For  this,  one  uses  the  theory  of  ply-elasticity 
and  the  methods  available  in  fracture  mechanics. 
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Accordingly,  Che  development  of  the  fracture  model  for  sub-laminate 
cracks  rests  upon  the  calculation  of  G(o,a)  and  the  physical  measurement 


of  G  . 
c 

The  Calculation  of  G(o,a) 

In  the  theory  of  the  LEFM,  the  singular  stress  field  near  a  crack  tip 
in  a  homogeneous,  isotropic  elastic  body  is  represented  by  analytical 
functions  in  the  theory  of  complex  variables  [26].  The  near-field  stresses 
are  obtained  for  three  particular  modes  of  cratck  extension.  These  are 
known  as  the  opening  mode  (I) ,  the  sliding  mode  (II)  and  the  anti-plane 
shearing  mode  (III).  For  each  mode,  the  stresses  are  expressed  in  terms 
of  the  associated  stress  intensity  factor  K  [27];  and  consequently,  the 
associated  G(o,a)  is  computed  in  terms  of  K.  Since  the  relation  between 
K  and  G  is  one-to-one.  Equation  (3)  reduces  to  the  form 

K(o,a)  -  Kc 

Similar  relations  between  K  and  G  for  orthotropic  media  having  a  crack 
orientated  along  one  of  the  major  axes  can  also  be  obtained  [28].  But  the 
analytical  solutions  for  the  singular  stress  field  often  require  tedious 
mathematical  derivations. 

Direct  solution  methods  for  G(o,a)  have  been  available;  among  them 
are  the  well-known  J-integral  method  [29,  30]  and  the  method  of  virtual 
crack  closure  technique  by  Irwin  [31]. 

Irwin  [31]  observed  that  the  elastic  strain  energy  released  during  a 
virtual  crack  extension  da  is  equal  to  the  work  done  in  closing  it  again. 
The  closing  of  Aa  provides  the  solution  for  the  surface  tractions  a 


over  Aa.  The  crack-tip  energy  release  rate  is  then  represented  by 

0  •  ikJT  <5-a=>d*  (5> 

where  Au  is  the  crack  opening  displacement  vector  over  Aa. 

If  the  crack  extension  involves  all  three  modes  (1,  II,  III),  the 
vector  product  in  Equation  (5)  will  give  a  sum  of  three  scalars,  associated 
respectively  with  G^,  G^,  and  ®III’ 

The  virtual  crack-closure  representation  is  particularly  adaptive  to 
machine  computations.  Rybicki  and  Kanninen  [32]  suggested  a  2-dimensional 
finite  element  technique  to  evaluate  G  for  a  line  crack  in  a  plane.  The 
crack-tip  stress  vector  cr  and  the  displacement  vector  Au  in  Equation  (5) 
are  approximated  by  the  nodal  forces  and  displacements  respectively,  in  a 
finite  element  representation  (for  detail,  see  [18]). 

Wang  and  Crossman  applied  this  technique  in  a  generalized  plane  strain 
finite  element  routine  [33],  which  can  simllate  a  line  crack  propagation 
in  the  2-dimensional  cross-section  of  a  laminate.  Since  a  general  lami¬ 
nate  under  load  may  suffer  cross-sectional  warping,  the  routine  actually 
computes  3-dimensional  stresses  and  displacements  [33]. 

If  a  laminate  is  subjected  to  the  far-field  stress  oq  and  a  sub¬ 
laminate  crack  is  to  be  simulated,  the  crack-tip  energy  release  rate  can 
be  expressed  in  the  general  form 

G  *  C  (a/t)'t*(a  /E  )2  (6) 

e  e  oo 

where  E  is  laminate  stiffness  in  o  direction,  and  t  is  a  characteristic 
o  o 

length,  taken  usually  as  the  thickness  of  the  ply.  Cg  is  a  function 

dependent  only  on  the  crack  size  a,  which  is  routinely  generated  for  a 
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given  type  of  crack  in  a  given  laminate. 

Similarly,  if  the  laminate  is  subjected  to  a  uniform  temperature  load 
AT,  and  if  a  thermally  induced  crack  opening  results ;  then  the  associa¬ 
ted  energy  release  rate  at  the  crack  tip  can  he  expressed  by, 


gt  “  cT(a/t)'t*(AT)i 


(7) 


Generally,  AT  is  characterized  to  represent  the  temperature  load  due  to 
laminate  curing;  but  it  can  include  subsequent  environmental e.g.  humidity) 
relaxation  effects  before  mechanical  load  is  applied  \_2Z],  The  combined 
temperature  and  mechanical  load  effects  on  the  stress  field  are  linearly 
additive;  and  the  total  available  energy  release  rate  is  expressed  as, 


G-  t(Cj1/2*ert  +  (C_)1/2-AT]2*t 
e  oi 


(8) 


where  e  is  the  far-field  strain  (■  cr  /E  ). 
o  o  o 

It  is  seen  that  the  finite  element  technique  is  extremely  versatile, 
and  can  be  efficiently  executed  for  simulating  complicated  sub-laminate 
cracking  such  as  delamination. 

Nonetheless,  the  accuracy  of  the  technique  has  been  a  subject  of  con¬ 
cern  of  many  analysts,  because  of  its  approximate  nature  in  representing 
a  mathematically  singular  stress  field.  As  has  been  demonstrated  by 
Raju  and  Crews  [34],  Spilker  and  Chou  [35],  and  Wang  and  Choi  [36],  the 
finite  element  stress  solutions  could  lose  significant  accuracy  in  the 
small  vicinity  of  the  singular  point;  the  region  in  which  the  stress  be¬ 
comes  Inaccurate  is  generally  much  smaller  than,  say,  a  fiber  diameter, 
due  to  the  very  nature  of  the  singularity  [36].  However,  it  can  also 
be  demonstrated  that  the  stresses  in  that  small  region  do  not  contribute 
significantly  to  the  crack  opening  energy  release  rate,  especially  for  a 
crack  size  much  larger  in  proportion. 
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Figure  2  shows  a  close  comparison  between  the  finite  element  computed 
G  and  the  elasticity  solution  counter-part,  for  a  transverse  crack  located 
in  the  mid-layer  of  a  3- layer  laminate  (see  inset  of  Figure  2).  The 
extreme  layers  are  designated  as  material  1  and  the  mid-layer  as  material 
2.  Both  were  assumed  homogeneous  and  Isotropic  materials.  The  exact  solu¬ 
tion  for  G  as  a  function  of  the  crack  size  a  is  given  by  Islda  [37],  using 
complex  stress  potentials;  while  the  finite  element  solution  is  computed 
using  a  rather  coarse  constant-strain,  triangular  element  mesh.  It  is 
seen  that  the  finite  element  solution  for  G(a)  compares  well  with  the  exact 
solution.  To  obtain  more  accurate  stresses  near  the  crack  tip, however,  one 
may  need  to  deploy  a  finer  element  mesh,  a  singular  element  or  a  higher- 
order  element  formulation  [38].  But,  such  extreme  measures  are  often  un¬ 
necessary  for  the  computation  of  G. 

The  Evaluation  of  G 

c 

When  fracture  occurs  in  the  material,  the  energy  required  in  the  pro¬ 
cess  is  expected  to  depend  on  the  morphology  of  the  fracture  surface, 
which  must  be  examined  at  the  microscale.  Fibrous  composites  are  known  to 
possess  complex  fracture  surface  details,  even  in  matrix-predominant  cracks 
The  observed  delamination  surfaces  in  graphite-epoxy  composites  show,  for 
Instance,  a  considerable  raggedness  because  the  crack  must  pass  around  the 
reinforcing  fibers.  But,  at  the  macroscopic  level,  the  crack  surface 
details  are  "smoothed  out";  and  their  effects  are  reflected  in  the  measured 
quantity  of  the  fracture  resistance  Gc> 

For  this  reason,  G^  measured  for  some  matrix-predominant  fractures  in 
epoxy-based  composites  has  been  found  to  depend  on  the  direction  of  frac¬ 
ture  propagation.  Cullen  [39]  and  Williams  [40]  considered  two  different 


cases  of  delamination  as  Illustrated  In  Figure  3.  The  first  case  is  0*/0® 

delamination  In  which  the  crack  path  Is  in  the  fiber  direction,  while  the 

second  is  90®/90°  delaalnation  where  the  crack  path  is  transverse  to  the 

fiber  direction.  The  experiments  were  performed  using  a  graphite-epoxy 

unidirectional  laainate,  subjected  to  mode-I  cracking  condition.  They 

found  that  the  aicroacopic  morphology  of  the  90®/90®  delaalnation  surface 

exhibited  considerably  more  raggedness  than  the  0®/0®  delaalnation  surface. 

This  resulted  in  marked  differences  for  the  measured  G  .  Note  that  these 

c 

two  fracture  events  occur  essentially  in  the  same  interface  when  viewed 
macro8copically.  Yet,  the  respective  values  can  differ  greatly  depend¬ 
ing  on  the  direction  of  the  crack  propagation. 

Hence,  when  the  Griffith  formula  (3)  is  applied  for  cracks  in  compo¬ 
sites,  the  term  G^  requires  a  precise  qualification.  Similarly,  when  a 

test  method  is  devised  to  measure  G  ,  it  is  also  necessary  to  consider  the 

c 

dimensional  and  directional  characteristics  of  the  measured  data. 

Mode-I  interlaminar  G^.  A  commonly  used  test  method  to  determine 
interlaminar  GIc  is  the  splitting  cantlliver  beam.  Cullen  [39]  and  Wilkins 
[41]  have  used  this  method  to  determine  the  Interlaminar  G_  when  the  crack 

IC 

is  propagating  in  the  direction  of  the  fibers  (0®/0®  delamination).  For 

some  graphite-epoxy  systems,  they  found  that  GJc  at  room  temperature  is 

2 

about  0.85  lb/ln,  or  130  J/m  .  This  value  is  about  twice  the  G,  measured 

ic 

for  pure  epoxy  resin.  Williams  [40]  used  a  compact  specimen  which  simulates 
roughly  a  90°/90®  delamination.  He  found,  for  the  same  material  system, 
a  GIc  value  of  1.3  lb/ in,  about  three  times  that  of  the  pure  resin.  Williams 
explained  that  the  fracture  surface  in  his  specimens  showed  fiber  breakage 
as  well  as  fractured  epoxy  debris;  this  had  resulted  in  a  higher  value 
for  Gjc  than  that  found  for  0®/0®  delamination  by  Cullen  [39]. 
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In  another  paper  by  Wilkins,  et.  al.  [42],  it  is  reported  that  G^c 
in  delamination  of  0°/90°  interface  is  also  higher  than  that  of  0°/0" 
delamination.  All  these  findings  indicate  the  directional  dependent  na¬ 
ture  of  Gjc 

Mixed-mode  interlaminar  G^  jjjc*  Th®  splitting  cantiliver  beam 
method  has  also  been  used  in  mixed-mode  cracking  experiment.  In  this  case, 
it  is  necessary  to  apply  different  loads  at  the  top  and  the  bottom  parts 
of  the  split  so  as  to  create  both  an  opening  (I)  and  a  sliding  (II) 
action.  Vanderkley  [43]  and  Wilkins  [41]  conducted  tests  on  the  same 
graphite-epoxy  system  (used  for  their  mode-I  tests) ,  and  found  the  total 
energy  release  rate  G^  n)c  "  Gic  +  Gnc  exhibited  a  strong  de¬ 

pendence  on  the  Gjj/Gj  ratio. 

This  phenomenon  is  not  uncommon  in  mixed-mode  fracture.  Similar 
observations  were  reported  for  brittle  metals  as  well  as  pure  epoxy  resins 
[44].  It  is,  perhaps,  more  pronounced  in  fibrous  composites.  Generally,  it 
IS  thought  that  the  increased  fracture  resistance  is  the  result  of  excess¬ 
ive  matrix  yielding  under  shear,  as  well  as  crack-closure  friction  due  to 
the  sliding  action. 

In  an  experiment  on  double-notched  off-axis  unidirectional  graphite- 

epoxy  laminates,  Wang,  et.  al.  [45]  measured  the  mixed-mode  G,T  as  a 

U-,  ii;c 

continuous  function  of  G__/G_,  see  Figure  4.  It  is  seen  that  G._  is 

II  I  (I,  II)c 

mono tonically  increasing  with  G^/G^.. 

But  in  the  same  experiment,  it  shows  also  that  the  opening  part  G 

ic 

remains  roughly  independent  of  G^/G^.  This  suggests  that  mixed-mode  crack 
is  essentially  controlled  by  mode-I.  Of  course,  such  a  suggestion  is  only 
temporary;  more  study  is  needed  to  fully  understand  the  true  nature  of  the 
mixed-mode  crack  mechanisms. 
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MECHANICS  OF  TRANSVERSE  CRACKS 


Physical  Mechanisms  at  the  Macroscale 

Transverse  cracks  are  found  In  epoxy-based  laminates,  even  at  a  low 

loading  level.  When  viewed  at  the  macroscopic  scale,  the  cracking  action 

is  simply  a  sudden  separation  of  fibers  by  breaking  the  epoxy  bond.  To 

illustrate,  consider  as  an  example  a  [0/90]  type  laminate  such  as  shown 

s 

by  the  inset  in  Figure  2.  There,  material  1  is  the  0°-layer  and  material 
2  is  the  90®-layer.  Under  the  far-field  tensile  load,  material  2  could 
suffer  multiple  transverse  cracks.  Generally,  the  sequence  of  events  is 
as  follows:  A  single  crack  forms  first  when  the  applied  load  reaches  a 
certain  critical  value,  which  defines  the  "onset"  of  the  events;  as  the 
applied  load  increases,  more  similar  cracks  are  formed.  If  there  is  no 
other  failure  mode  setting  in  at  high  load  (e.g.  0°-rupture,  delamination, 
etc.},  the  number  of  transverse  cracks  will  continue  increasing,  until  it 
reaches  a  saturation  density. 

Figure  5  shows  a  load-sequence  x-radiographs  taken  for  a  [0/90]g 
graphite-epoxy  laminate  under  ascending  tensile  load.  Transverse  cracks 
in  the  90°-layer  are  seen  to  form  in  increasing  numbers.  From  these  photo¬ 
graphs,  a  plot  of  crack  density  (e.g.  cracks  per  inch  of  specimen  length) 
versus  the  applied  load  can  be  obtained.  Figure  6  shows  a  family  of  such 

experimental  plots  for  a  series  of  [0/90  /0]  laminates,  n  ■  1,  2,  3,  4. 

n 

Examination  of  the  plots  in  Figure  6  reveals  several  interesting  fea¬ 
tures.  First,  the  onset  load  for  the  first  transverse  crack  seems  to  be 
Influenced  profoundly  by  the  90°-layer  thickness,  especially  when  it  is 
very  thin.  Take,  for  example,  the  case  of  n  *  1  in  Figure  6,  the  onset 
load  is  almost  twice  that  of  n  ■  2. 


Secondly,  the  crack  density  also  shows  dependence  on  the  same  thickness 
factor.  Generally,  the  laminate  with  thinner  90#-layer  Is  capable  of  yield¬ 
ing  a  higher  crack  density.  (But,  for  the  case  of  n  ■  1,  failure  of  0*-layer 
at  high  load  Interrupted  the  development  of  more  transverse  cracks ) . 

The  90°-layer  thickness  effect  on  transverse  cracking  was  first  docu¬ 
mented  experimentally  by  Bader,  et.  al.  [46].  They  attributed  the  effect 
to  the  constraining  actions  of  the  adjacent  0°-layers.  Observing  that  a 
transverse  crack  can  be  no  larger  than  the  90°-layer  thickness,  the  energy 
release  rate  which  drives  the  crack  is  thus  limited  by  the  same  factor. 

It  is  the  total  strain  energy  trapped  in  the  90°-layer  which  determines  the 
onset  of  the  cracking,  not  the  ln-situ  tensile  stress. 

As  for  the  effect  on  crack  density,  it  has  been  explained  by  the  exist¬ 
ence  of  a  "shear-lag"  zone  at  the  transverse  crack  root  [46  -  48] .  That  is 
an  interlaminar  shear  stress  is  developed  on  the  0/90  interface  where  a 
transverse  crack  terminates.  This  shear  stress  is  singular  at  the  crack 
root,  but  decays  exponentially  a  distance  away  [20].  Similarly,  the  in-sltu 
tensile  stress  in  the  90°-layer  is  nil  at  the  transverse  crack,  but  it  re¬ 
gains  its  far-field  value  outside  the  shear-lag  zone.  Thus,  ideally,  any 
two  adjacent  cracks  should  be  spaced  by  the  shear-lag  distance.  Since 
this  distance  is  proportional  to  the  90°-layer  thickness  [20],  hence  the 
observed  thickness  effect  on  crack  density  [46]. 

Although  the  shear-lag  concept  is  ideally  correct,  the  so-called 
"characteristic"  spacing  of  transverse  cracks  does  not  occur  exactly  in  prac¬ 
tice,  due  to  the  random  nature  of  material  flaws.  However,  transverse  cracks 
may  cause  other  failure  modes,  such  as  delamination  C2l].  Failure  modes  inter¬ 
actions  will  be  discussed  later  in  this  paper. 
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AC  Che  microscopic  level,  Che  mechanisms  of  a  cransverse  crack  are 
much  more  perplexing.  For  example,  Che  exacC  kinemacics  of  crack  formaclon 
Is  noc  fully  known.  Post-test  SEM  examination  of  Che  crack  surfaces  gen¬ 
erally  gives  a  very  ragged  appearance;  ciny  epoxy  debris  and  sometimes 
broken  fibers  are  also  seen  [40] .  Figure  7  shows  an  x-ray  plane  view  of  a 

[+25/90„]  laminate  after  cransverse  cracking,  left,  and  an  edge-view 
—  is 

micrograph  on  the  right.  It  is  seen  chat  the  cransverse  crack  is  practi¬ 
cally  a  plane  crack  of  a  rectangular  dimension,  which  is  bounded  by  Che 
width  of  the  specimen  and  the  90°-layer  thickness.  There  is  no  evidence 
Co  indicate  chat  Che  crack  formation  is  progressive  in  nature.  In  fact, 
all  experiments  tend  to  suggest  a  sudden  dynamical  formation. 

This  dynamical  nature.  Chough  still  conjectured,  has  been  qualified 
by  many  who  monitored  the  acoustics  emitted  during  the  crack  formation 
(see,  e.g.  [49]). 

While  It  is  difficult  to  reduce  these  physically  observed  facts  into 
a  general  law,  they  nevertheless  provide  the  necessary  rationale  for  the 
formation  of  an  analytical  model.  In  what  follows,  an  energy  formulation 
is  presented,  which  describes  the  most  essential  observed  characteristics 
of  the  transverse  cracking  phenomenon. 

The  Energy  Formulation 

For  purpose  of  clarity,  consider  a  [0/90]  type  laminate,  as  shown  in 

s 

Figure  8.  It  will  be  assumed  that,  in  the  90s-layer,  the  material  has  a 
random  distribution  of  microflaws.  The  gross  effects  of  the  microflaws  at 
the  macroscopic  scale  are  represented  by  a  characteristic  distribution  of 
"effective"  flaws  which  cannot  be  physically  seen.  But,  under  stress, 
these  effective  flaws  are  capable  of  propagating  suddenly  into  transverse 
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cracks,  which  are  physically  real.  The  individual  size  of  the  effective 
flaws  is  denoted  by  2a;  and  any  two  adjacent  flaws  are  spaced  by  a  dis¬ 
tance  S,  see  Figure  8. 

Given  a  unit  length  of  the  specimen,  there  is  a  probability  density 
function  f(a),  and  a  probability  density  function  f(S).  For  sake  of  no 
evidence  to  suggest  otherwise,  the  two  functions  are  assumed  to  take  a 
form  for  normal  distributions  [50]; 

£(»)  ■  exp  [-(a  -  vJ2/2o2]  (9) 

*< <s>  *  s7!?  <-<s  -  “a’2'20,2)  <10> 

where  p  and  o  are  the  mean  and  the  standard  deviation  of  the  respective 
distribution  functions. 

Among  the  effective  flaws,  the  size  of  the  "worst"  one  is  denoted  by 

2a  .  For  definiteness,  assume  a  the  99  percentile  of  f(a).  That  is  99Z 
o  o 

of  the  flaws  are  smaller  than  2a  . 

o 

Then,  under  the  far-field  load,  say  e  ,  the  "worst"  flaw  2a  will 

o  o 

become  the  first  transverse  crack.  The  critical  value  for  e  at  the  onset 

o 

is  calculated  by  substituting  Equation  (8)  into  the  Griffith  criterion  (3) 
[(Ce)1/2-  eo  +  (CT>1/2.  AT]2-  t  -  Gc  (11) 

where  C  and  C_  are  evaluated  at  a  ■  a  . 
e  T  o 

As  has  been  detailed  earlier,  the  energy  release  rate  coefficients 
C£  and  C^,  are  generated  numerically  by  the  finite  element  crack-closure 
routine,  given  the  geometry  and  material  moduli  of  the  laminas.  Thus, 
in  Equation  (11),  all  quantities  except  e  will  be  given:  AT  is  the  tern- 


perature  load  due  to  cooling  and  possibly  also  environmental  relaxation;  t 

is  a  characteristic  length( thickness  of  one  ply)  and  G  Is  measured  experl- 

c 

mentally  simulating  mode-I  90°/90°  delamlnation(see  Figure  3-b) . 

Now,  the  difficulty  rests  upon  the  choice  of  a  ,  or  for  the  same  matter 

o 

the  choice  of  p  and  a  in  the  distribution  functions  in  (9). 
a  a 

Clearly,  2aQ  must  be  smaller  than  the  908-layer. thickness.  But  the 
latter  can  be  made  arbitrarily  large.  Hence,  a  finite  bound  on  a0  exists 
even  if  the  90°-layer  thickness  is  unbounded. 

As  an  example,  consider  the  experimental  results  reported  in  [21]. 

The  tensile  strength  of  T300/934  [90g]  laminates  averaged  -  7000  psi. 
Assume  that  the  "largest"  effective  flaw  in  the  £9003  laminate  has  a  size 
2aQ(paralell  to  the  fiber)}  this  flaw  will  determine  the  laminate's  tensile 
strength  given  by  the  Griffith  formula  (2).  From  which  one  obtains 

a  -  G  E/ir a  2  (12) 

o  c  u 

With  E  ■  1.7  X  10^  psl  and  G  ■  0.9  lb/in  for  0°/0°  mode-I  delamination, 

c 

Equation  (12)  determines  a  -  0.01",  or  about  2  times  the  ply  thickness  of 

o 

the  T300/934  systems. 

It  may  thus  be  inferred  that  the  largest  flaw  in  an  unbounded  90° 

laminate  (or  In  a  £o/90  3  type  laminate  with  n  *  2)  should  be  no  greater 

n  s 

than  ao~0.01".  Indeed,  in  several  experimental  correlations  conducted 

by  Wang  and  Crossman  [19,  20,  21],  using  the  same  material  system,  a  is 

o 

found  in  the  order  of  0.0075",  or  1.5  times  the  ply  thickness. 

Numerical  Examples  by  Monte-Carlo  Simulation 


considered  laminates  are  in  the  form  [0./90  ]  ,  n  «  1,  2  and  3.  The 

l  ns 

material  system  is  the  AS-3501-06  graphite-epoxy  system  (for  the  material 
moduli,  see  [51]).  The  nominal  ply  thickness  of  this  system  is  0.0052". 

In  order  to  define  the  parameters  y  and  a  in  the  f(a)  and  f(S)  distri 
bution  functions,  the  following  values  are  chosen  for  the  laminates  n  *  2 
and  3: 

y  -  0.0036",  o  -  0.0013"  (13) 

a  a 


y  -  0.0125",  0  -  0.0046"  (14) 

s  s 

These  are  chosen  to  reflect  the  fact  that  a  ,  being  over  99%  of  f(a), 

o 

may  take  a  value  (see  [50]), 

a  -  y  +  3o  _  0.0075"  (15) 

o  a  a 

and  that  for  y  *  0.0125"  it  implies  80  effective  flaws  to  the  inch.  The 
s 

choice  of  the  standard  deviations  is  a  matter  of  adjusting  to  the  actual 
data  scatter  [51]. 

As  for  the  laminates  of  n  ■  1,  the  thickness  of  the  90°-layer  is  only 

0.0104",  or  a  <  0.0052".  For  this  case,  the  choice  of  y  and  o  are  as 
o  a  a 

follows: 

y  -  0.0021",  o  -  0.0C075"  (16) 

a  a 


so  that 


y  +  30  ,  0.0044" 


(17) 


The  choice  of  the  spacing  parameters  remains  the  same  for  n  •  1 


The  Monte-Carlo  simulation  procedures  start  with  the  generation  of  a 
set  of  N  (  «  80)  flaws  whose  random  sizes  and  spaclngs  are  represented  by 
the  respective  density  functions  (9)  and  (10).  This  is  done  by  generating 
first  a  set  of  N  random  values  in  the  interval  (0,1).  Then,  by  equating 
the  cumulative  function  of  f(a)  to  each  of  the  random  values,  a  random 
set  {a^},  i  *  1,  N  is  computed.  Among  the  values  in  {a^},  the  probability 
of  the  largest  to  be  equal  or  greater  than  aQ  as  defined  in  (15)  or  (17) 
is  about  99%. 

Similarly,  a  random  set  {S^},  i  ■  1,  (N  -  1) ,  is  also  generated  by 
the  random  number  scheme.  is  then  assigned  to  be  the  spacing  between 
the  ith  and  the  (i  +  1)  th  flaws. 


Thus,  a  computer  research  follows,  which  determines  the  flaw  most 
likely  to  become  a  transverse  crack.  The  first  to  occur,  clearly,  is  the 
worst  flaw  in  the  {a^}  set,  and  the  corresponding  applied  load  for  the 
onset  of  the  first  crack  is  then  determined  using  Equation  (11). 

After  the  first  crack  is  formed,  another  flaw  in  the  {a^}  set  will  be 
come  a  crack  at  a  slightly  higher  load.  But  this  flaw  is  located  at  a 
random  distance  from  the  first  crack.  And  the  presence  of  the  physical 
(real)  crack  has  a  stress  reducing  effect  on  the  rest  of  the  flaws.  Thus, 
if  the  second  flaw  to  become  a  crack  is  located  outside  the  shear-lag 
zone,  it  will  not  feel  the  presence  of  the  first  crack;  and  hence,  its 
available  energy  release  rate  at  the  instant  of  cracking  is  given  by 
GQ(a,a),  the  same  as  in  Equation  (8). 

On  the  other  hand,  if  the  second  flaw  is  located  inside  the  shear-lag 


zone,  then  the  available  energy  release  rate  is  reduced  by  a  factor  de- 


pending  on  Its  distance  from  the  first  crack. 

Generally,  for  the  flaw  to  become  the  kth  crack,  one  must  search  to 
the  left  and  to  the  right  for  the  nearest  cracks.  If  so,  the  energy  release 
rate  at  the  instance  of  the  kth  crack  is  given  by 

Gk(a,a)  -  R(S£)  GQ(a,a)  R(Sr>  (18) 

where  S.  is  the  distance  to  the  left  crack,  S  the  distance  to  the  right 

A*  IT 

crack. 

In  Equation  (18),  R(S)  is  the  energy  rate  retention  factor  which  is 

generated  by  the  finite  element  routine  for  a  flaw  of  2a  (<2a  )  and  is 

o 

placed  at  varying  distance  S  from  a  transverse  crack,  see  Figure  9. 

For  the  laminate  family  (0„/90  ]  ,  a  single  R(S)  curve  can  be  genera- 

l  ns 

ted  if  expressed  in  terms  of  S/nt,  t  being  the  thickness  of  one  90°-ply, 
see  Figure  9.  Note  that  100%  energy  rate  retention  is  expected  beyond 
S  ■  9nt;  the  latter  is  actually  the  size  of  the  shear-lag  zone  [19]. 

The  applied  load,  corresponding  to  the  kth  crack  is  determined  from 

Gfc(a,a)  -  Gc  (19) 

The  search  essentially  simulates  the  cracking  process  as  it  would 
occur  naturally.  Each  random  simulation  represents  an  actual  test  case; 
and  repeated  simulations  represent  actual  tests  on  replicas.  Table  A 
illustrates  the  flow-chart  for  the  Monte-Carlo  simulation. 

Figures  10,  11,  and  12  show  the  simulated  crack-density  versus  load 
plots  for  respectively,  [02/90]s,  [<?2/ 902 3  and  [O2/9O3]  laminates.  In  the 
simulations,  the  quantities  appearing  in  Equation  (11)  are  assigned  as 


follows : 


(2 


AT  -  225°F,  t  -  0.0052",  G  -  1.3  lb/in 

c 


The  shaded  band  in  each  of  these  plots  is  the  corresponding  experi 
mental  data  band  from  4  test  specimens.  These  test  data  were  reported 
earlier  in  Reference  [52];  and  the  details  for  the  simulation  computer 
routine  are  reported  in  Reference  [51]. 

MECHANICS  OF  FREE  EDGE  DELAMINATION 
The  Classical  Free-Edge  Problem 

The  free  edge  delamination  problem  has  attracted  increased  research 
interest  ever  since  the  advent  of  composite  laminates.  The  reason  is  that 
delamination  is  one  of  the  most  damaging  sub-laminate  failure  modes.  Gen¬ 
erally,  it  is  a  plane  crack  which  forms  along  the  laminate  free  edge  and 
propagates  inward  along  an  interface  of  two  adjacent  layers. 

'  Figure  13  shows  an  x-ray  plane  view  of  a  [+45/0/90]  graphite-epoxy 

3 

laminate  under  uniaxial  tension  (left).  Delamination  is  seen  to  occur 
along  both  edges  of  the  laminate,  with  essentially  uniform  growth  toward 
the  center.  An  edge  view  of  this  crack  is  shown  on  the  right,  which  indi¬ 
cated  that  the  crack  is  basically  contained  inside  the  90°-layer,  not 
necessarily  in  any  one  given  layer  interface.  Moreover,  the  cracked  plane 
is  quite  zig-zagged  along  the  length  of  the  free  edge. 

Similar  photographs,  taken  for  a  [+25/90.  ,„]  laminate  after  failure, 

—  1/  l  s 

are  shown  in  Figure  14.  Here,  the  edge  crack  is  seen  to  have  formed  in¬ 
side  the  90°-layer,  but  branched  to  the  25/90  interface  as  it  propagated. 
The  branch-out  is  due  to  the  skewed  cracks  which  occur  in  the  90°-layer 
ahead  of  the  edge  de lamination  front. 


Note  that  both  laminates  have  a  90°-layer;  and  they  could  suffer  trans¬ 
verse  cracking  under  the  tensile  load.  But,  because  the  thickness  of  the 
90°-layer  is  thin,  a  high  load  is  required  for  transverse  cracking  accord¬ 
ing  to  the  energy  theory  presented  earlier.  Instead,  in  this  case,  edge 
delamination  is  induced  as  the  first  event  of  sub-laminate  failure. 

Analytical  studies  of  the  free  edge  delamination  problem  have  origina¬ 
ted  from  the  much  celebrated  work  of  Pipes  and  Pagano  [53],  who  formulated 
and  computed  the  boundary-layer  interlaminar  stress  solutions  for  a  long 
symmetric  laminate  under  tension.  The  formulation  is  based  on  the  macro¬ 
scopic  ply-elasticity  theory,  which  regards  the  material  layers  as  indi¬ 
vidually  homogeneous  media;  material  and  geometrical  discontinuities  exist 
only  across  the  layer  Interfaces. 

Generally,  the  free  edge  stress  field  is  three-dimensional,  and  is 
singular  at  the  intersection  of  the  free  edge  with  a  layer  interface  Q36]. 
Hence,  interface  delamination  is  caused  by  the  highly  concentrated  edge 
stresses,  especially  the  interlaminar  normal  stress  o^. 

In  a  series  of  tensile  strength  tests,  Bjeletich,  et.  al.  [54] 

examined  the  failure  modes  of  six  families  of  quasl-lsotropic  laminates  by 

alternating  the  stacking  sequence  of  the  0°,  90°  and  +45°  layers.  Edge 

stress  analysis  indicated  a  compressive  a  along  the  free  edge  of  the 

[0/90/+45]  laminate,  while  a  tensile  a  for  the  [+45/0/90]  laminate.  The 
—  s  z  —  s 

latter  developed  extensive  delamination;  and  the  growth  of  it  had  led  to 
a  much  lowered  tensile  strength.  Clearly,  knowledge  of  the  free-edge 
stresses  can  provide  an  explanation  why  delamination  occur;  but  a  quanti¬ 
tative  prediction  for  its  occurrence  requires  a  more  precisely  defined 


criterion. 


Ia  a  study  on  delamination  for  similar  graphite-epoxy  laminates, 

Rodini,  et.  al.  [55]  determined  experimentally  the  critical  tensile  loads 

at  the  onset  of  free  edge  delamination  in  a  [+45  / 0/90  ]  ,  n  ■  1,  2,  3 

—  n  n  ns 

family.  They  found  that  the  critical  laminate  tensile  stress  o^  varied 

greatly  with  the  value  of  n.  Specifically,  the  critical  o^  decreases  at 

the  rate  of  about  ViT,  even  though  an  edge  stress  analysis  yields  identical 

o  for  the  same  a  for  all  values  of  n. 
z  x 

Such  a  90°-layer  thickness  dependent  behavior  is  related  to  that  found 
in  the  transverse  cracking  problems,  suggesting  the  observed  phenomenon 
is  again  fracture  in  nature. 

The  Energy  Criterion 

In  the  work  of  Wang  and  Crossman  [18],  it  is  assumed  that  material 
flaws  exist  randomly  on  any  one  of  the  interfaces  between  the  material 
layers.  Those  flaws  which  are  located  within  the  free-edge  stress  zone 
form  an  "effective"  flaw  having  a  size  aQ  at  the  instance  of  onset  of  de- 
lamination,  see  inset  in  Figure  15.  Following  the  finite  element  crack- 
closure  procedure  by  allowing  virtual  extension  of  aQ  along  the  layer  in¬ 
terface,  an  energy  release  rate  curve  is  generated,  such  as  shown  in 
Figure  15. 

It  is  seen  that  the  energy  release  rate  G  increases  sharply  with  a, 

but  reaches  an  asymptote  at  a«  a  .  Generally,  the  value  of  a  is  about 

m  m 

one-half  the  layer-thickness  which  contains  the  delamination  [18].  Thus, 
the  layer  thickness  affect  the  value  of  a  ,  and  also  the  value  of  G. 

IQ 

.Expressing  the  computed  energy  release  rate.  Equation  (6) , 


G(e  ,  a)  -  C  VC  ( a)  AT]2  t 


(21) 


and  applying  the  Griffith  criterion  (3)  for  the  onset  of  delaaination,  one 
obtains  the  critical  far-field  strain, 

<«0>cr  "  C  'fcj*  *  AT]  /  ,  a  =  a0  (22) 

Note  that  in  (22),  the  thermal  residual  effect  is  also  included.  But 
such  thermal  effect  on  free  edge  delaaination  is  generally  minimal  [~56]. 

The  problem  comes  down  to  two  important  questions;  namely,  what  value 
is  aQ,  and  which  Interface  is  likely  to  delaminate? 

To  answer  the  first  question,  recall  the  earlier  discussions  about  a^ 
in  unbounded  90°  laminates.  As  inferred  from  the  Griffith's  criterion,  a 

o 

is  in  the  order  of  0.01",  or  2  times  the  ply  thickness  for  commercial 
graphite-epoxy  systems.  It  is  believed  that  aQ  for  delamination  is  at 
least  of  this  magnitude  if  not  larger,  because  of  possible  additional  cut¬ 
ting  flaws  along  the  free  edge.  In  any  event,  if  a^  in  Figure  15  is  less 
than  2  times  of  the  ply  thickness,  one  simply  uses  a  ■  a^  and  predicts 
from  Equation  (22)  the  minimum  possible  load  for  the  onset  of  delamina¬ 
tion.  Thus 

(eo>cr  “  E  ^7*  *  ^7*^  .  a  =  affi  (23) 

The  answer  to  the  second  question,  however,  is  complicated,  if  not 
lncondse.  Consider  as  an  example  a  under  uniform  compres¬ 

sion.  Interlaminar  tensile  a is  developed  in  this  laminate  and  edge  de- 
lamination  is  Induced  before  global  buckling,  see  Figure  16.  A  3-dimen¬ 
sional  through-thickness  display  of  a  is  shown  in  Figure  17(a).  It  is  seen 

z 

that  the  largest  a  is  located  on  the  45/ —45  interface  (in  fact,  this  a 


is  singular  at  the  free  edge).  On  the  other  hand,  a  along  the  mid-plane 

z 

(z  ■  0)  is  finite  at  the  free  edge.  In  addition,  there  is  also  a  singular 

interlaminar  shear  stress  t  along  the  45/-45  interface.  Thus,  the  stress 

xz 

analysis  suggests  the  45/-4S  interface  as  the  most  probable  delamination  site. 

However,  according  to  the  energy  release  rate  argument  as  well  as 
the  "effective"  flaw  hypothesis,  one  must  calculate  the  energy  release  rate 
curves  along  all  possible  Interfaces  in  order  to  determine  the  interface 
most  energetically  favorable  for  delamination. 

The  finite  element  crack-closure  results  of  the  G(a)  curves  for  this 
example  problem  is  shown  in  Figure  17(b).  Actually,  for 'a  >  t,  the  mid¬ 
plane  interface  yields  a  much  larger  energy  release  rate  than  the  45/-4 5 
interface.  For  a  <  t  however,  the  energy  release  rate  on  45/-45  interface  is 
slightly  larger  than  that  on  the  mid-plane.  For  all  likelihood,  the  worst 
interface  flaw  size  aQ  is  larger  than  tj  and  the  most  energetically  favor¬ 
able  delamination  site  should  be  the  mid-plane. 

This  conclusion  is  confirmed  by  experiments  conducted  by  Wang  and 
Slomiana  [57],  who  also  investigated  compression-induced  delamination  in 
[0_/90_/+45,]  and  [0/90/0/ 90/+45/+45]  laminates.  In  all  cases,  mid- 
plane  and  mode-I  delamination  was  predicted  and  observed.  Using  the  energy 

criterion  of  Equation  (23),  and  setting  G  ■  1.0  lb/in,  the  predicted 

c 

minimum  onset  loads  (c^)  for  these  laminates  are  tabulated  below,  along 
with  the  corresponding  test  da ta( average  of  3  specimens  each) > 


In  the  previous  examples,  edge  delamination  is  the  first  and  only  sub¬ 
laminate  failure  mode  before  laminate  buckling.  The  onset  load  for  delami¬ 
nation  is  synonymous  to  the  final  laminate  failure  load.  In  most  other 
cases,  edge  delamination  and  transverse  cracking  are  often  two  competing 
failure  modes,  which  interact  each  other  through  a  complicated  mechanisms. 

Examine  again  the  photomicrographs  shown  in  Figures  13  and  14,  where 
a  [+45/0/90)s  and  a  [+25/90^2^8  l*®*11®1®8  were  lo*ded  by  axial  tension. 

The  thickness  of  the  90°-layer  in  both  laminates  is  so  small  that  trans¬ 
verse  cracking  is  possible  only  at  a  high  load.  Instead,  an  edge  delami¬ 
nation  is  induced  as  the  first  sub-laminate  failure  mode.  Nonetheless,  the 
stress  field  in  the  90°-layer,  which  contains  the  delamination,  appears  to 
be  extremely  complex.  The  skewed  cracks  shown  in  Figure  14  are  indicative 
of  the  complex  stress  field  ahead  of  the  delamlnetion  front.  But,  just 
how  profoundly  these  secondary  cracking  events  influence  the  delamination 
growth  is  far  from  known. 

On  the  other  hand.  Figure  18(a)  shows  an  x-ray  plane  view  of  a  [+25/902) 
laminate  after  delamination.  Actually,  this  laminate  suffered  first  trans¬ 
verse  cracks  in  the  90°-layer;  subsequent  edge  delamination  in  the  25/90 
Interface  developed  at  a  higher  load.  It  is  seen  from  the  picture  that 
delamination  actually  caused  many  more  transverse  cracks  because  of  stress 
concentration  along  the  curved  delamination  front.  On  the  other  hand,  in  the 
undelaminated  protion  only  a  few  transverse  cracks  are  in  fact  present. 

This  phenomenon  has  repeatedly  been  observed  in  experiments  [21];  but  no 
serious  account  has  been  attempted  to  analytically  model  the  mechanisms. 


Law  [56]  assumed  that  when  transverse  cracks  occur  in  the  90°-layer 


the  tensile  modulus  E^  and  the  Poisson  ratio  of  the  90°-layer  both 

reduce  to  nearly  zero.  With  reduced  moduli  in  the  90°-layer,  the  overall 

laminate  becomes  energetically  delamination  prone  in  the  25/90  interface, 

resulting  in  a  mixed-mode  cracking  (I,  II).  Some  details  of  Law's  analysis 

will  be  discussed  later  in  a  case  study. 

Figure  18(b)  shown  an  x-ray  plane  view  for  a  [+25/904]a  laminate  soon 

after  transverse  cracking  took  place.  In  this  case,  scattered  transverse 

cracks  appeared  first,  which  immediately  caused  delamination  on  both  25/90 

interfaces  in  the  area  where  they  occurred.  The  delamination  is,  in  fact, 

a  combined  effect  stemming  from  the  transverse  crack-tip  and  free  edge 

stress  concentrations.  Figure  19  illustrates  this  combined  effect  in  an 

isometric  view.  The  interface  shear  stress  x  ,  which  is  developed  at  the 

xz 

transverse  crack  root,  tends  to  delaminate  the  25/90  interface  in  the  x- 

direction;  and  the  interlaminar  normal  stress,  o  ,  which  is  a  free  edge 

z 

effect,  tends  to  open  up  the  same  interface  in  the  y-dlrectlon.  The  com¬ 
bined  action  gives  rise  to  a  mixed-mode,  delta-shaped  delamination.  Evi¬ 
dence  of  this  type  of  delamination  can  be  found  by  examining  Figure  18b. 

Law  [56]  approximated  the  combined  effect  by  first  computing  for 
the  x-direction  delamination  by  a  2-dimensional  analysis,  and  then  by 
computing  for  the  y-direction  as  in  the  case  of  free  edge  delamlnatlon. 

Since  both  G  and  G  are  considered  crack-driving  forces,  their  vector 
x  y 

sum  is  used  to  predict  the  onset  of  the  delta-shaped  delamination.  Ob¬ 
viously,  such  an  approximate  approach  is  only  tentative;  the  problem 
remains  to  be  treated  more  rigorously  by  a  truly  3-<iiniensional  analysis 
simulating  a  "contoured"  plane  delamination.  But,  as  it  will  be  shown  in 
the  next  section,  Law's  approach  was  at  least  qualitatively  correct. 
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It  is  now  apparent  that  sub- laminate  cracking  ,  especially  delamination, 
are  physically  complicated  even  under  simple  tension.  The  energy  method  just 
presented  should  be  applicable  if  the  interaction  mechanisms  amongst  different 
sub-laminate  crack  modes  could  be  more  accurately  simulated,  this  topic,  how¬ 
ever,  remains  to  be  further .developed. 

An  Experimental  Case  Study 

An  experimental  case  study  on  transverse  cracking,  free  edge  delamina¬ 
tion  and  their  mutual  Interaction  was  conducted  by  Crossman  and  Wang  [21]. 

They  chose  a  family  of  graphite-epoxy  laminates  (T300/934)  in  the  form  of 

[+25/90  ]  .  By  varying  the  single  parameter  n  -  1/2,  1,  2,  3,  4,  6  and  8, 

—  n  s 

they  were  able  to  control  the  occurrence  of  the  various  competing  failure 
modes  in  the  laminates  when  loaded  in  tension.  Both  x-radiography  and 
microphotography  were  used  to  monitor  the  sub-laminate  cracking  development 
as  a  function  of  the  applied  load.  For  clarity,  a  graphical  summary  of 
their  experimental  findings  is  shown  in  Figure  20.  For  n  *  1/2  and  1,  the 
laminates  are  classified  as  having  a  "thin"  90°-layer.  Transverse  crack¬ 
ing  is  suppressed  due  to  the  lack  of  sufficient  strain  energy  stored  in 
the  90°-layer.  Hence,  free  edge  delamination  occurs  as  the  first  sub¬ 
laminate  failure  mode.  For  n  •  1/2,  a  stable  mode-I  delamination  is  pro¬ 
duced  during  the  loading  until  final  failure.  For  n  ■  1,  however,  trans¬ 
verse  cracks  are  induced  in  the  area  of  delamination.  Because  of  these 
cracks,  the  delamination  plane  is  often  branched  to  the  25/90  Interface. 

This  latter  appearance  is  indicative  of  an  interaction  taking  place  between 
the  two  failure  modes. 

When  n  »  2  and  3,  the  laminates  are  denoted  as  having  "thick"  90°- 
layer.  In  both,  transverse  cracks  occur  first  long  before  delamination. 

It  may  be  thought  that  the  onset  of  transverse  cracking  is  an  independent 
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event,  while  delamination  is  caused,  at  least  partially,  by  the  presence 
o£  the  multiple  transverse  cracks.  Indeed,  the  delamination  plane  in  the 
case  of  n  ■  2  resembles  that  for  n  •  1;  and  the  delaminatlon  plane  in  the 
n  ■  3  laminate  is  entirely  in  the  25/90  interface.  Also,  for  n  *  3,  the 
growth  of  delamination  becomes  more  rapid. 

The  laminates  with  "very  thick"  90#-layer  are  those  for  n  ■  4,  6, 
and  8.  In  all  cases  transverse  cracks  occur  first  followed  immediately 
by  delta-shaped  delaminations  emanating  from  the  transverse  crack  roots. 
These  soon  coalesce  to  form  large  scale  25/90  delamination,  resulting  in 
premature  laminate  failure. 

Figure  21(a)  summarises  the  various  onset  loads  Sx  for  each  of  the 
laminates  (signified  by  the  value  of  n).  Since  the  laminates'  longitu¬ 
dinal  Young  modulus  varies  with  the  value  of  n,  It  Is  sometimes  more 

convenient  to  use  the  laminate  tensile  strain,  e  as  a  measure  of  the 

x 

applied  load.  Note  that  except  for  n  *  1/2  and  1,  transverse  cracking  is 
the  first  sub-laminate  failure  event.  Generally,  growth  of  multiple  cracks 
is  fully  developed,  followed  by  interlayer  delamination.  However,  as 
n  >  3,  the  load-gap  between  the  transverse  cracking  and  delamination  events 
grows  closer.  In  fact  for  n  ■  8,  these  two  events  are  closely  linked  to 
the  final  laminate  failure. 

Note  also  that  except  for  n  ■  1/2,  1  and  2,  onset  of  delaminatlon  and 
final  failure  are  separated  only  by  less  than  5%  of  loading. 

By  means  of  the  energy  model  presented  earlier.  Law  [56]  conducted  a 
computer  simulation  of  the  transverse  cracking  and  delamination  processes 
in  the  (+25/90nls  family.  The  numerical  results  are  shown  in  Figure  21(b), 
which  is  superscribed  on  the  experimental  results  shown  in  Figure  21(a). 


Law  has  actually  conducted  four-  cypes  of  crack  simulations:  (a)  onset 
of  transverse  cracking  in  the  90°-layers  for  all  values  of  n;  (b)  mid-plane 
(z  ■  0),  mode-I  delamination  for  n  *  1/2,  1,  2,  and  3;  (c)  25/90  interface 
delamination  for  n  ■  2,  3,  and  4  (the  90°-layer  is  saturated  with  trans¬ 
verse  cracks;  and  is  represented  by  reduced  moduli);  and  (d)  transverse 
crack  root/ free  edge  combined  delamination  of  the  25/90  interface  for  n  ■  3, 
4,  6  and  8  (as  illustrated  in  Figure  19). 

It  is  seen  from  Figure  21b  that  the  predicted  onset  loads  for  transverse 
cracking  agree  well  with  the  experiment,  while  the  predicted  delamination 
loads  (represented  by  3  curves)  form  an  envelope  to  which  the  experimental 
points  fall  below.  Although  the  predicted  delamination  loads  for  n  ■  6  and 
8  are  off  considerably  from  the  experimental  points,  the  energy  analysis 
seems  to  capture  the  physical  mechanisms  of  these  various  failure  modes; 
the  analysis  gives  correctly  the  trend  of  the  failure  loads. 

Final  Failure  Mechanisms 

In  the  case  study  discussed  above,  the  sub-laminate  failure  sequence 

is  seen  to  lead  to  the  final  rupture  of  the  laminate.  It  is  generally 

thought  that  final  failure  is  determined  essentially  by  the  strength  of 

the  load-carrying  plies;  i.e.  in  the  case  of  the  [+25/90  ]  family,  the 

n  s 

+25-plies  should  carry  practically  all  the  load.  But,  what  determines  the 
strength  of  the  +25-plies?  And,  is  it  adversely  influenced  by  the  sub¬ 
laminate  cracking  events  that  proceded  the  final  laminate  rupture?  These 
questions  relate  to  the  final  failure  modes  and  their  associated  fracture 
mechanisms. 

Crossman  and  Wang  C583  conducted  an  extensive  experimental  investiga¬ 
tion  into  the  failure  sequence  of  several  families  of  laminates  under  uni- 


axial  tension.  Some  of  his  results  may  help  to  answer  the  above  questions. 

Consider  first  the  effect  of  thickness  of  the  load-carrying  plies  by 
examining  the  results  in  Table  B. 

It  is  noted  that  the  final  failure  mode  of  [+25/+2S]  is  due  simply 

~  1  "  3 

to  fiber  breaking  across  the  laminate,  while  the  final  failure  mode  of 

[25_/-25_]  involves  a  major  transverse  crack  in  the  center  core  (-25.°- 
2  Z  S  4 

layer)  followed  by  multiple  fiber-splitting  in  the  252~layer.  In  particu¬ 
lar,  there  is  no  fiber  breakage  involved  in  the  latter  case,  see  Figure  22. 

An  energy  analysis  of  the  respective  failure  modes  indicated  that  the 
transverse  cracking  mode  in  the  [+25/+251  is  suppressed  by  the  actual 

~  S 

thickness  of  the  25°-layer;  there  is  not  enough  strain  energy  to  cause  this 
failure  mode  at  a  lower  load  level.  Hence,  it  forces  a  fiber  breakage 
instead.  However,  by  doubling  the  layer  thickness,  in  the  case  of  [25^1-25^] 
the  transverse  cracking  mode  emerges  at  a  much  lower  load.  The  difference 
in  the  failure  modes  results  in  the  difference  in  the  applied  ultimate  loads. 

Clearly,  the  strength  of  the  so-called  "load-carrying"  plies  in  a  lami¬ 
nate  depends  also  on  the  particular  fracture  modes  that  cause  failure.  In 
the  examples  above,  the  ply-thickness  factor  is  seen  again  to  play  a 
profound  role  in  the  final  failure  mechanisms. 

When  the  respective  load-carrying  plies  are  laminated  with  a  90°-layer, 
such  as  in  [+25/+25/902 ] s  and  [ 252/" 252/ 90^ ] g »  transverse  cracks  in  the 
90°-layer  and  edge  delamination  are  found  to  proceed  any  failure  in  the 
25°-plies.  From  the  results  displayed  in  Table  B,  it  is  seen  that  the 
sub-laminate  cracking  events  neither  changed  the  failure  mode  of  the 
j^S'-plies,  nor  decreased  the  final  ultimate  loading,  although  the 
presence  of  the  90°-layer  in  [ 252/-252/ 902 ] g  slightly  increased  the  lami- 


nate  load-carrying  capacity.  The  reason  is  that  [252/-252]g  has  four  (4) 

25®-plies  in  the  core,  while  [25^1  a  ^as  on^y  two  25°-plies 

scattered  by  the  presence  of  the  90°-layer  in  the  core. 

However,  adverse  effects  of  sub-laminate  cracking  on  laminate  final 

failure  modes  are  seen  from  the  test  results  displayed  in  Table  C.  Here, 

the  failure  modes  in  [+0/90- ]  and  [+0/90-1  are  compared.  Note  that  the 

—  is  —  os 

90°-layer  thickness  in  the  latter  is  exaggerated  in  order  to  separate 
widely  the  respective  sub-laminate  crack  modes. 

In  the  cases  of  [0-/90-]  and  [+25/90- ]  ,  the  final  failure  mode  in 

l  0  s  —  OS 

the  load-carrying  plies  involves  multiple  fiber-splitting  in  0°  or  25°  layers 
over  the  root  of  a  90°-layer  transverse  crack;  see  Figure  22.  On  the 
other  hand,  there  is  no  fiber-splitting  in  the  °r  t*ie  Ci^S/ 90j  1  g 

laminates. 

A  three  dimensional  stress  analysis  by  Vang,  et.al.Q59l  has  shown  that 
over  the  90°-layer  transverse  crack  root,  stress  concentration  can  be  amp¬ 
lified  several  times  by  the  presence  of  the  laminate  free  edge.  In 

particular,  the  size  of  the  stress  concentration  zone  is  proportional  to 
the  90°-layer  thickness;  which  explains  why  the  considerable  difference 

in  the  final  loads  between  [0o/90„]  and  [0„/90o]  ;  and  between  [+25/ 90_ ] 

i  Z  s  Z  0  s  —  l  s 

and  [+25/90-]  ,  see  Table  C. 

—  os 

Fiber-splitting  of  the  outer  ply  seems  to  trigger  the  final  rupture  of  the 

[0o/90o]  and  [+25/ 90_ ]  laminates.  The  splitting  mechanisms  is  due  to 
zos  —  os 

the  exlstance  of  the  tensile  stress  normal  to  the  fibers,  which  is  magni¬ 
fied  in  amplitude  by  the  90®-layer  transverse  crack  root. 

Similar  stress  condition  also  exists  in  the  (+30/90- ]  laminates.  But 

—  os 

in  this  case,  the  tensile  stress  (normal  to  +30°-fibers)  is  not  large  enough 


to  cause  fiber  splitting;  see  also  Table  C.  As  a  result,  the  final  failure 
inodes  for  [+30/90,,]  and  [+30/90.]  all  involve  fiber  breaking;  and  their 
final  failure  load  are  essentially  the  same.  Thus,  a  mere  increase  in  the 
angle  from  25°  to  30°  changes  the  failure  modes  in  C±25/90q!18  and.  C*30/90q38« 
Clearly,  the  final  failure  modes  in  the  "load-carrying"  plies  can  be 
influenced  by  nature  from  the  sub-laminate  cracking  .  Moreover,  there  are 
probably  a  host  of  other  possible  modes  that  are  competing  for  dominance. 

The  specific  failure  occurrence  in  a  given  laminate  is  generally  determined 
by  the  lamination  geometry  as  well  as  the  intrinsic  material  properties. 

And,  of  course,  there  are  always  the  statistical  uncertainties  in  the  "in¬ 
trinsic"  material  properties. 


CONCLUSIONS 

This  paper  has  attempted  to  give  an  overview  of  a  fracture  mechanics 
approach  to  some  of  the  frequently  observed  sub-laminate  cracks  in  epoxy- 
based  composite  laminates.  The  laminates  discussed  in  the  paper  served 
mainly  to  illustrate  how  a  certain  failure  sequence  is  developed 
and  how  the  predictive  model  is  constructed.  They  are  certainly  not  the  ones 
which  may  be  used  in  practical  design  applications.  The  major  objective  here 
is  to  enable  to  see  through  the  various  fracture  mechanisms  and  try  to 
understand  them;  to  this  end,  one  often  had  to  exaggerate  the  geometric  para¬ 
meters  (e.g.  the  thickness  of  90°-layer)  in  order  to  ascertain  their  full  range 
of  influence. 

Indeed,  the  field  of  fracture  in  composites  seems  infinitely  complex. 

Yet,  numerous  past  and  present  investigators  have  toiled  the  field  with 
many  successes.  And,  new  discoveries  are  being  continuously  made  to 
further  its  advancement.  The  results  reported  in  this  paper  represent 
but  one  effort  to  fill  the  dimensional  span  from  micromechanics  to  struc¬ 
tural  mechanics  in  composite  failure  analysis. 
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Mi crograph( right)  Showing  Two  Transverse  Cracks  In  A  C±2 5/90  ,,3  Laminate 
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Energy  Release  Rate  Retention  Factor  R(S)  for  a  Flaw  Located  at 
Distance  S  from  a  Transverse  Crack. 
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Comparison  of  Experiment  and  the  Monte-Carlo  simulations  of  Transverse  Cracks 
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Comparison  of  Experiment  and  the  Monte-Carlo  Simulations  of  Transverse  Cracks 
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Cob pari son  of  Experiment  and  the  Monte-Carlo  Simulations  of  Transverse  Cracks 
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Figure  17  Interlaminar  Normal  Stress  Distribution  and  Energv  Felease 

Rate  Coefficients  for  Under  Compression.  Ref.C57l 


T300/934  (+25/90n)$,  n=t/2,l ,2, 3,4, 6, 8. 


<a>  (b)  .  (c) 


SCHEMATIC  OF  THE  FRACTURE  SEQUENCE  IN  THE  (+25/90  )  LAMINATES 

—  ns 

(a)  JUST  PRIOR  TO  EDGE  DELAMINATION,  (b)  SUBSEQUENT  TO  EDGE 
DELAMINATION,  (c)  JUST  PRIOR  to  FINAL  FAILURE.  Ref.[2l]. 
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Figure.  22  Final  Failure  Modes:  (a)  Cross-sectional  Fiber 

Break;  (b)  25-layer  Transverse  Cracks;  (c)  0-layer 
Fiber  Splits  over  90-layer  Transverse  Cracks. 
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ABSTRACT 

This  paper  examines  the  sequence  of  sub  laminate 

fracture  and  the  damage  state  which  exists  just  prior  to 

laminate  failure  under  tensile  loading.  Tests  were  conducted 

on  T300/934  symmetric  laminates  of  [theta, -theta,  90 1  type, 

s 

in  which  the  thickness  of  the  90  degree  ply  varies  from  2  to 
16  plies  and  the  angle  "theta"  takes  on  values  of  0,  25,  30, 
and  60  degrees.  Incremental  tensile  loading  followed  by  DIB 
enhanced  x-radiography  is  the  primary  means  of  documentation. 
Deplying  experiments  and  real-time  x-ray  radiography  during 
tensile  loading  are  also  described.  Failure  sequence  and  mode 
are  shown  to  change  with  stacking  sequence  and  ply  thickness 
in  a  manner  which  supports  a  model  of  fracture  based  on  a 
through-the-thickness  local  redistribution  of  stress  to  the 
primary  load  carrying  ply. 


INTRODUCTION 


This  paper  presents  a  detailed  investigation  of  the  sub-laminate 
fracture  modes  which  exist  just  prior  to  laminate  fracture  in  a  series  of 
laminates  of  the  [theta/-theta/90]g  family  where  the  angle  "theta"  takes 
on  values  of  0,  25,  30,  60  degrees.  This  paper  has  been  preceded  by  a  series 
of  experimental  and  analytical  studies  which  have  addressed  the  fracture 
of  unnotched  tensile  coupons  of  the  T300/934  graphite-epoxy  material.  Pre¬ 
vious  experimental  investigations  have  included  the  documentation  of 
fracture  mode  and  sequence  of  fracture  in  quasi-isotropic  [Q/90/45/-45] 

s 

laminates  by  metallographic  means  [1]  and  the  description  of  fracture 

sequence  in  [25/-25/90]  laminates  by  radio  opaque  dye  penetration  of  edge 

s 

cracks  and  x-radiography  [2]. 

Transverse  cracking  of  the  90  deg  and  off  axis  plies  and  delamination 
at  the  free  edges  of  the  coupon  were  observed  to  occur  at  tensile  loads 
sufficiently  lower  than  the  ultimate  strength  that  these  fracture  modes  could 
be  viewed  as  important  precursors  to  laminate  fracture.  Analytical  methods 
based  on  strain  energy  release  rate  [3]  and  stochastic  processes  [5]  have 
been  successful  in  predicting  the  onset  and  accumulation  of  transverse  cracks 
under  both  static  and  cyclic  loading.  The  utility  of  strain  energy  release 
rate  analysis  to  predict  the  onset  of  free-edge  delamination  and  local  de¬ 
lamination  induced  by  transverse  cracking  has  also  been  established  under 
tensile  loading  conditions  [4],  However,  more  recent  analyses  of  delami¬ 
nation  under  fatigue  cycling  [6,7]  have  demonstrated  that  delamination  growth 
under  mixed  mode  crack  opening  is  not  simply  related  to  knowledge  of  delami¬ 
nation  initiation  and  growth  under  Mode  I  crack  opening  and  static  loading 


conditions. 


The  detailed  analysis  of  transverse  cracking  and  isolated  delaminations 
provide  examples  of  what  Reifsnider  and  colleagues  have  termed  the  "inductive" 
approach  to  composite  failure  analysis  [8].  This  approach  requires  that  the 
analytical  models  account  for  the  detailed  experimentally  observed  micro- 
structural  fracture  processes  for  their  interaction.  A  detailed  review  and 
historical  perspective  of  the  experimental  observations  and  analytical  modeling 
of  laminate  fracture  tinder  static  and  cyclic  loads  is  provided  in  [8]. 

An  important  corollary  to  the  inductive  approach  is  that  SOME  of  the 
observed  sub-laminate  fracture  modes  may  not  have  to  be  modeled  and  analyzed 
in  detail  if  the  prediction  of  final  laminate  fracture  is  of  primary  concern 
[12],  What  is  needed  is  a  fracture  model  which  accounts  for  the  development 
of  the  sub-laminate  fracture  condition  ("damage  state")  which  exists  just 
prior  to  fracture  and  an  understanding  of  why  that  condition  is  critical. 

The  purpose  of  this  paper  is  to  document  the  sequence  of  sub-laminate 
fracture  processes  and  to  describe  the  final  damage  state  of  a  series  of 
simple  laminate  constructions.  The  results  of  this  study  indicate  a  strong 
dependence  of  that  final  state  on  fiber  angle,  ply  thickness,  and  stacking 
sequence.  With  few  exceptions,  the  final  fracture  event  is  triggered  not  by 
a  sharp  stress  concentration  but  by  a  stress  redistribution  between  plies  in 
a  region  that  is  small  compared  to  the  size  of  the  coupon. 

EXPERIMENTAL  PROCEDURE 

Laminates  were  fabricated  from  Fiberite  T300/934  prepreg  tape  and  auto¬ 
clave  cured  at  450  K  according  to  the  manufacturer's  recommended  cure  proce¬ 
dure.  Nominal  fiber  volume  ranged  from  64  to  68  percent,  and  ply  thickness 
ranged  from  .0119  to  .0127  cm.  The  30  cm  by  30  cm  panels  were  cut  by 
diamond  saw  to  a  coupon  length  of  22.5  cm  and  a  width  of  2.5  cm.  Fiber¬ 
glass  tabs  of  length  3.75  cm  were  bonded  to  each  end  of  the  coupon.  Coupons 


were  scored  ac  ambient  laboratory  conditions  (20-23  deg.  C  and  55-65  percent 
humidity)  and  tested  6-12  months  after  fabrication.  All  tests  were  conducted 
on  an  Instron  Test  Machine  at  a  controlled  displacement  rate  of  25mm/mm.  An 
average  ultimate  tensile  strength  for  each  type  of  laminate  was  determined  by 
ramp  tension  testing  of  two  replicas.  Stress  strain  curves  were  based  on  the 
measured  output  of  1/4  in.  EA-13-125TM-120  type  strain  gages. 

Three  to  five  replicas  of  each  laminate  type  were  step  loaded  to  approxi¬ 
mately  50,  70,  80,  and  90  percent  of  the  average  UTS,  and  then,  based  on  the 
audible  cracking  which  the  experienced  operator  observed,  the  coupon  was  step 
loaded  in  progressively  smaller  increments  to  failure.  The  first  load  step  was 
chosen  to  coincide  with  the  onset  of  transverse  cracking  as  determined  in  pre¬ 
vious  studies  [2,  11].  At  least  two  coupons  of  each  series  were  step  loaded  to 
between  97-99  percent  of  UTS  and  examined  prior  to  final  failure.  The  total 
number  of  step  loadings  varied  between  5  and  10  for  a  given  coupon. 

The  examination  after  each  step  loading  and  unloading  consisted  of  (1)  re¬ 
moving  the  coupon  from  the  tensile  tester,  (2)  swabbing  the  edges  of  the  coupon 
with  di-iodobutane  (DIB)  solution  to  penetrate  the  cracks,  (3)  wiping  off  the 
excess  DIB  after  30  s,  and  (4)  placing  the  coupon  atop  a  Polaroid  film  (Type  55 
Positive/Negative  ASA  50)  in  a  Faxatron  804  point  source  x-ray  chamber  for  a  1.25 
min.  exposure  at  30  KV.  Following  the  x-ray  exposure  the  specimen  was  optically 
examined  on  edge  at  30X  to  document  the  through-the-thickness  location  of  ply 
cracks  and  delamination. 

Several  coupons  of  [0/0/90  ]  (n*l,  2,  4,  8)  and  [25/-25/90  ]  (n*8)  were 

ii  s  ns 

deplled  following  tensile  loading  to  70  percent  of  UTS.  The  deplying  procedure 
developed  by  Freeman  [9]  was  used  and  consisted  of  heating  to  700  K  for  1.25 
hrs  to  facilitate  the  separation  of  the  plies.  Deplied  0  and  25  deg.  layers 
were  then  examined  under  scanning  electron  microscopy  (SEM)  for  damage  to  the 


fibers  in  Che  load  carrying  plies. 

Several  coupons  of  [25/-25/90  ]  (n-8)  type  were  examined  under  real-time 

n  s 

x-ray  exposure  during  tensile  loading.  The  Lockheed  HERTIS  (High  Energy 
Real  Time  Inspection  System)  facility  was  used  for  this  purpose.  The  system 
consisted  of  a  45KV/10mA  soft  x-ray  source  with  a  4mm  focal  point  on  one  side 
of  the  coupon  being  loaded  in  a  portable  5000  lb  (22.  2kN)  hydraulic  test 
frame  and  fluorescent  screen  and  vidicon  tube  to  pick  up  the  x-ray  image  on 
the  other  side  of  the  coupon.  At  increments  of  20  percent  of  the  specimen  UTS, 
the  specimen  was  unloaded  and  the  coupon  edges  were  coated  with  DIB  X-ray  en¬ 
hancing  penetrant.  Two  vidicon  tubes  capable  of  imaging  a  1  by  1  cm  and  12  by 
12  cm  area  respectively  were  used  in  the  experiments.  The  vidicon  signal  was 
then  digitally  processed  and  signal  enhanced  for  viewing  in  real  time  and  on 
a  video  tape  system. 

EXPERIMENTAL  RESULTS 

Method  of  Presentation 

The  cataloging  of  x-ray  and  optical  observations  of  sub-laminate  fracture 
of  the  laminate  series  tested  in  this  program  resulted  in  hundreds  of  individual 
photographs  which  provide  a  snapshot  of  damage  accumulation  and  which,  taken  as 
a  whole,  provide  some  measure  of  failure  mode  and  sequence  variability  within 
a  given  laminate  construction. 

In  the  presentation  of  these  results  we  have  chosen  to  represent  the 
damage  state  obtained  just  prior  to  laminate  failure  by  a  set  of  schematic 
drawings  and  graphs  which  succinctly  compare  the  accumulation  of  damage  within 
a  given  series  of  similar  laminates. 

Fig.  1  presents  a  summary  of  fracture  sequence  in  [0/0/90n]g  (n-1,  2,  4,  8) 
laminates.  The  diagram  labeled  (A)  in  Fig.  1(a)  represents  the  center  portion 
of  the  gage  length  of  the  n«l  type  tensile  coupon  as  recorded  on  the  DIB  enhanced 


x-radiograph.  The  tensile  axis  is  vertical  and  the  width  direction  is  hori¬ 
zontal  in  this  view. 

The  diagram  labeled  (B)  in  Fig.  1(a)  depicts  a  view  of  the  specimen 
edge.  The  tensile  axis  is  vertical  and  the  through-the-thickness  direction 
is  horizontal  in  this  view.  The  dotted  lines  indicate  ply  boundaries  at  the 
(0/0)  and  (0/90)  ply  interfaces.  In  subsequent  schematics  the  (theta/-theta) 
and  (-theta/ 90)  boundaries  are  dotted. 

The  third  diagram  labeled  (C)  in  Fig.  1(a)  presents  an  optical  view  of 
the  gage  length  of  the  specimen  after  laminate  fracture.  The  tensile  axis  is 
vertical  and  the  width  direction  is  horizontal  in  this  view. 

To  the  right  of  these  three  diagrams  is  a  graph  in  histogram  form  con¬ 
taining  the  range  of  onset  stresses  for  each  of  several  individual  sub-laminate 
fracture  modes.  All  four  graphs  in  Fig.  1  have  the  same  horizontal  stress 
scale  to  facilitate  a  comparison  of  fracture  sequence  and  to  provide  the  de¬ 
veloper  of  fracture  process  models  with  some  sense  of  the  interaction  between 
fracture  modes  within  a  given  family  of  laminates.  The  onset  data  for  a  given 
fracture  nssde  is  displayed  at  progressively  higher  vertical  positions  on  each 
graph  according  to  the  order  in  which  each  occurs  during  tensile  loading. 

Co  /90  1  Laminates. 
n'  nr's _ 

Fig.  1  summarizes  the  sub-laminate  fracture  sequence  in  a  series  of  [0/0/ 

90  ]  laminates  where  n*l,2,4,8.  The  damage  noted  by  (1)  in  Fig.  1(a)  denotes 
transverse  cracking  which  occurs  early  during  tensile  loading  of  all  four  lam¬ 
inates.  The  damage  mode  noted  by  (2)  in  Fig.  1(b)  indicates  the  presence  of 
local  delamination  at  the  [0/90]  interface  which  accompanies  the  initiation  of 
90  deg  transverse  cracks.  These  local  delaminations  are  observed  to  initiate 
more  frequently  at  stresses  approaching  the  UTS  of  the  laminate.  The  first 
transverse  cracks  which  initiate  do  not  exhibit  this  characteristic  delamination 


The  damage  near  note  (3)  in  Fig.  1(d)  denotes  the  presence  of  transverse 

cracks  in  the  0  deg.  ply.  Under  tensile  loading  only  the  n*8  laminates  were 

found  to  exhibit  this  ply  splitting  mode.  Reifsnider  et  al.  [8]  have  found 

splitting  to  be  much  more  common  in  [0/90]  type  laminates  during  fatigue 

s 

cycling. 

Examination  of  the  final  fracture  diagrams  in  Fig.  1  indicates  that  the 
amount  of  0  deg.  layer  splintering  increases  as  n  increases.  This  observation 
is  consistant  with  the  increased  delamination  between  the  0  and  90  deg.  plies 
for  increasing  n.  The  n=8  shattering  of  the  0  deg.  plies  closely  resembles  the. 
fracture  mode  of  unidirectional  0  deg.  tensile  specimens.  With  the  exception 
of  the  0  deg.  splitting  in  the  n*8  laminate,  there  is  no  obvious  precursor  to 
failure  of  the  load  carrying  ply.  The  initiation  and  accumulation  of  trans¬ 
verse  cracks  in  the  90  deg.  ply  appears  to  have  little  influence  on  the  failure 
of  the  0  deg  layer  in  all  of  the  other  laminates. 

Fig.  2  presents  a  summary  of  incremental  tensile  tests  conducted  on  a 

general  series  of  [0  /90  ]  laminates  in  which  the  0  deg.  layer  was  located  at 

n  ms 

either  the  surface  (a,c,d)  or  midplane  (b) .  While  the  onset  for  transverse 
cracking  is  clearly  influenced  by  temperature  (due  to  altered  thermal  residual 
stresses  [11]),  the  results  in  Fig.  2(c)  show  no  similar  temperature  dependence 
of  UTS.  Indeed,  the  only  case  where  a  lower  than  expected  UTS  occurs  is  in  the 
testing  of  the  [0/90/0]  laminate.  As  shown  in  Fig.  2(d)  it  is  the  only  laminate 
in  which  transverse  cracking  is  not  observed  prior  to  failure. 

Fiber  damage  has  been  observed  in  deplied  0  deg.  layers  at  the  tip  of  trans¬ 
verse  90  deg.  cracks  in  fatigue  loaded  specimens  [8].  In  this  sutdy  [0/0/90n]s 
n-8  laminates  were  deplied  after  tensile  loading  to  70  percent  of  UTS  and  ex¬ 
amined  for  fiber  damage.  The  results  shown  in  the  Fig.  3  SEM  photographs  in¬ 
dicate  that  far  less  than  one  percent  of  0  deg.  f.'bers  are  broken  in  the  vicinity 


of  transverse  cracks  under  static  loading.  The  possiblity  exists,  however 


that  the  iniatiation  of  a  new  90  deg.  crack  at  higher  stresses  could  trigger 
the  sufficient  fiber  damage  that  the  0  deg.  layer  fails  prematurely,  since 
the  energy  release  rate  from  transverse  cracking  increases  with  the  square  of 
the  stress  [3]. 

(25/-25/901  ,n=2,8  Laminates 

_ _ n-'s  _ 

Fig.  4  summarizes  the  sub-laminate  fracture  of  [25/-25/90  ]  laminates 

ns 

tested  in  this  program.  A  more  extensive  summary  including  n»l/2,l,3,4,6  lami¬ 
nate  types  is  found  [2].  The  region  marked  by  (1)  in  Fig.  4(a)  denotes  the 
prevalent  free— edge  delamination  which  is  at  the  -25/90  interface.  Delamination 
initiates  at  load  very  close  to  UTS  but  incremental  loading  of  these  laminates 
shows  that  it  is  the  initiation  of  25  deg.  transverse  cracking  denoted  by  (2) 
in  Fig.  4(b)  that  triggers  the  final  failure. 

In  Fig.  5(a)  one  can  observe  dark  horizontal  bands  on  the  25  deg.  fiber 
layer  in  locations  where  x-rays  taken  prior  to  failure  indicated  the  presence 
of  transverse  cracks  and  local  -25/90  interface  delamination.  These  regions 
are  marked  by  (1)  in  Fig.  5(a).  Examination  of  the  25  deg.  ply  surface  in 
these  regions  by  SEM  shows  that  the  bands  are  caused  by  the  pattern  of  resin 
adherence  and  do  not  show  any  damage  to  the  fibers.  This  obserbation  was  ver¬ 
ified  by  examination  of  deplied  n=8  laminates  as  well. 

HERTIS  real  time  x-ray  experiments  were  conducted  on  [25/-25/90n]g  n«8 
laminates  loaded  in  tension.  In  Fig.  5(d)  a  video  Image  of  one  specimen  at  99 
percent  of  ultimate  load  indicates  the  presence  of  a  90  deg.  transverse  crack 
and  an  associated  [-25/90]  local  delamination  at  location  (2)  in  the  figure. 

Stop  action  playback  of  the  video  taped  failure  of  this  specimen  showed 
that  the  final  fracture  occured  at  this  location  by  the  propagation  of  a  crack 
transverse  to  the  tensile  axis  from  the  edge  toward  the  center  of  the  specimen. 
After  propagating  by  fiber  breakage  to  about  half  way  across  the  specimen  the 


remaining  25  deg.  layers  failed  by  shearing  parallel  to  the  25  deg.  fiber 
axis.  Thus  the  critical  event  appears  to  involve  sufficient  load  redistribu¬ 
tion  into  the  uncracked  25  deg.  ply  that  fiber  failure  can  be  triggered  and 
propagate  transverse  to  the  tensile  axis. 

[30/ -30/^0  ]  and  [60/-60/W  1  Laminates 
n  s  n  s  _ 


Figs.  6  and  7  summarize  the  results  of  incremental  tensile  testing  of 
theta*  30  deg.  and  theta*  60  deg.  laminates.  Transverse  cracking  in  these  lami¬ 
nates  occurs  irregularly  and  most  often  in  association  with  the  initiation  of 
another  fracture  mode  as  at  location  (1)  in  Figs.  6  and  7.  The  transverse 
cracking  of  the  30  degree  plies  near  the  free  edge  of  the  laminate  is  found  to 
accompany  the  initiation  of  free-edge  delamination  at  the  -30/90  interface 
as  at  locations  (2)  in  Fig.  6.  In  many  cases  the  transverse  30  degree  crack 
causes  a  local  delamination  at  the  30/-30  interface  and  then  transverse  crack¬ 
ing  of  the  -30  deg.  ply.  With  the  exception  of  90  deg.  transverse  cracking, 
all  of  these  events  occur  very  near  the  UTS  of  this  laminate.  The  highest 
density  of  90  degree  transverse  cracks  are  located  under  the  delaminated  edges. 

In  n=l  laminates  transverse  cracking  is  vertually  absent  prior  to  failure. 
Furthermore,  the  delamination  plane  runs  randomly  through  the  90  deg.  layer, 
while  for  n*2,4,8  the  -30/90  interface  is  the  preferred  plane. 

In  the  theta*30  deg.  laminates,  the  final  fracture  is  by  propagation  of  a 
fiber  breaking  crack  transverse  to  the  tensile  axis  accompanied  by  a  small  per¬ 
centage  of  shear  out  of  30  deg.  layers  near  the  edge  of  the  coupon. 

In  theta*60  deg.  laminates,  the  critical  fracture  event  changes  from  the 
initiation  of  a  60  deg.  transverse  crack  for  n*l,2  to  the  initiation  of  a  90  deg 
transverse  crack  in  n=4,8  laminates.  In  either  case  the  presence  of  the  first 

transverse  crack  in  one  ply  triggers  the  initiation  of  several  transverse 
cracks  in  the  other  type  of  ply. 

In  the  n«4  laminate  the  90  deg.  crack  triggers  only  internal  splitting 


of  the  -60  and  60  layers  which  are  observed  after  the  DIB  penetrates  to  these 
cracks  along  the  90  transverse  crack  which  is  open  to  the  coupon  edge.  This  is 
noted  at  location  (2)  in  Fig.  7. 

For  n*l,  2  the  final  failure  occurs  by  propagation  of  the  transverse  crack 
in  the  outer  60  deg.  ply  and  a  fiber  breaking  crack  transverse  to  the  tensile 
load  in  the  -60  deg.  ply.  In  n=4,  8  the  failure  is  almost  entirely  due  to  shear 
of  the  60  deg.  layers  with  very  little  fiber  breakage. 

Closer  examination  of  the  fracture  patterns  in  Figs.  4,  6,  and  7  shows  that 
the  distribution  of  transverse  cracking  and  delamination  along  the  gage  length 
of  the  specimen  is  very  irregular.  It  is  unlikely  that  strain  measurement  via 
a  single  small  strain  gage  or  global  strain  measurement  via  extensometers  span¬ 
ning  the  gage  length  can  provide  sufficient  information  regarding  the  local  re¬ 
distribution  of  stresses  due  to  the  interaction  of  these  fracture  modes.  The 

"stiffness"  changes  associated  with  the  [theta/-theta/9G]  family  of  laminates 

s 

are  clearly  highly  dependent  on  the  location  and  size  of  the  zone  over  which 
strain  is  measured. 

Effect  of  90  deg.  Ply  Thickness 

If  stiffness  measurements  provide  little  of  value  in  assessing  the  fracture 
criterion  for  these  [theta,  -theta,  90 3 g  laminates,  we  are  left  with  a  need  to 
detail  the  degree  of  local  damage  and  determine  whether  the  undamaged  plies  at  a 
given  cross-section  of  the  laminate  can  carry  the  entire  tensile  load.  It  is 
important  to  understand  that  cracked  plies  can  still  carry  some  load  even  in  a 
direction  normal  to  the  crack  plane.  The  density  of  transverse  cracks  in  the 
region  of  greatest  damage  in  the  [theta,  -theta,  90]g  laminates  is  high  just 
prior  to  laminate  fracture.  Fig.  8  plots  the  inverse  of  crack  density  normalized 
to  the  thickness  of  the  90  deg.  ply.  Most  of  the  laminates  tested  have  S/T 


values  between  1-2. 


Finite  element  analyses  [10,  12]  of  cracked  90  deg.  layers  in  a  [0/90]g 
laminate  have  shown  that  the  effective  transverse  modulus  of  the  90  deg.  layer 
is  reduced  to  0.65  of  its  initial  value  at  S/T  =  2.0  and  to  0.40  at  S/T  **  1.0. 
Furthermore,  if  the  transverse  cracks  cause  a  local  delamination  at  the  theta/ 
90  interface  the  90  deg.  ply  is  more  fully  isolated  from  the  theta  load  car¬ 
rying  plies. 

On  the  basis  of  these  results  and  the  observed  high  density  of  transverse 
90  cracks,  it  is  realistic  to  consider  that  the  tensile  load  carried  by  the  90 
deg.  plies  is  very  small  compared  to  that  carried  by  the  theta  plies.  This 
assumption  is  borne  out  when  the  ultimate  tensile  load  is  plotted  for  each  of 
the  laminate  sequences  in  this  study,  as  shown  in  Fig.  9.  For  n*l,  2,  4  there 
is  very  little  dependence  of  ultimate  load  on  the  thickness  of  the  90  deg.  ply. 
For  theta=  0,  25  deg.  and  n*8  there  is  a  drop  in  load  carrying  ability,  while 
for  theta=  30,  60  and  n=8,  the  load  carrying  ability  is  slightly  enhanced.  The 
existance  of  large  regions  of  locally  delannated  theta/90  interface  in  the 
theta=  0,  25  deg.  and  n=8  laminates  is  a  possible  cause  of  the  observed  degra¬ 
dation. 


To  further  investigate  the  concept  of  local  stress  redistribution  as  a 

model  for  predicting  final  failure,  a  series  of  laminates  were  prepared  from 

the  [25/-25/90]  family  and  incrementally  tensile  tested  by  the  procedures 
s' 

already  described.  The  laminate  sequences  tested  are  shown  in  Fig.  10(a-d). 
Because  of  differences  in  stacking  sequence,  the  edge  view  diagrams  for  cases 
(a,  c,  d)  must  be  explained.  In  Fig.  1(a)  the  layup  has  no  90  deg.  ply  and 
the  transverse  cracks  drawn  across  the  center  layer  in  the  diagram  are  simply 
to  connect  the  transverse  cracks  which  are  continuous  across  the  inner  -25 
layers.  In  Fig.  10(c,  d)  the  outer  layer  is  taken  to  be  only  the  outer  25  dec 
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ply,  the  next  dotted  layer  actually  contalna  the  next  three  25  type  plies,  and 
the  center  layer  is  the  90  ply. 

The  layups  in  Fig.  10(a)  and  10(b)  contain  double  layers  of  25  deg.  plies 
adjacent  to  each  other.  In  an  earlier  study  [2]  the  doubling  up  of  25  deg. 
plies  resulted  in  altering  the  failure  node  from  one  of  fiber  breakage  to  one 
of  transverse  cracking  and  shear  out  parallel  to  the  25  deg.  fiber  axis. 

These  results  were  duplicated  here  for  cases  (a)  and  (b) . 

In  the  [25/25/— 25/— 25Js  laminate,  the  critical  fracture  event  was  the 
initiation  of  a  single  transverse  crack  in  the  outer  25  deg.  ply.  This  was 
followed  immediately  by  a  fine  transverse  splitting  of  the  inner  two  plies  of 
-25  deg.  and  subsequent  shear  out  of  the  remaining  25  deg.  ply. 

The  [25/25/-25/-25/90/90]g  laminate  developed  a  high  density  of  transverse 
90  deg.  ply  cracks  early  in  the  tensile  loading,  followed  by  the  initiation  of 
-25/90  delaminatlon  at  the  free  edge  and  a  rapid  unstable  propagation  of  de- 
lamlnatlon  toward  the  specimen  centerline  from  both  edges.  The  speed  with 
which  the  delaminatlon  growth  takes  place  is  evident  in  the  wavlness  of  the 
final  delaminatlon  front  near  the  centerline  of  the  specimen.  Earlier  studies 
[3,  4]  have  shown  that  the  driving  force  (strain  energy  release  rate)  drops  off 
rapidly  as  the  delaminatlon  approaches  the  center  line  thus  allowing  this  region 
to  remain  undelaminated  until  final  laminate  fracture. 

The  stacking  sequences  in  (c)  and  (d)  show  the  same  early  failure  modes 
of  transverse  cracking  of  the  90  deg.  ply  and  25/90  interface  delamination. 

As  in  laminate  (b)  the  outer  25  degree  ply  tends  to  transverse  crack  near  the 
coupon  free  edges.  However,  in  the  case  of  (c)  and  (d)  a  local  delamination 
is  formed  between  the  outer  25  deg.  ply  and  the  inner  25  deg.  plies  and  trans¬ 
verse  crack  do  not  form  in  the  adjacent  -25  deg.  ply.  Because  the  cracks  in 
the  outer  25  deg.  ply  are  isolated  from  the  adjacent  25  deg.  plies  by  local 


delamination,  further  ply  degradation  is  inhibited  until  another  5-10  KSI  is 
carried  by  the  laminate.  As  a  result,  laminates  (c)  and  (d)  have  a  higher  UTS 
than  laminate  (b)  and  exhibit  a  mode  of  failure  dominated  by  fiber  breakage 
rather  than  shear  out. 

Fig.  11  provides  a  summary  of  the  range  of  ultimate  loads  carried  by  the 
four  laminates  described  in  Fig.  10.  If  one  compares  the  [25/25/-25/-2S] 

s 

results  for  [25/25/-25/-25/90/90]  ,  the  addition  of  the  90  layer  contributes 

s 

little  to  the  tensile  load  carrying  ability  of  the  laminate,  as  observed  in 

the  [theta,  -theta,  90]  series.  However,  when  the  25/-25  plies  are  alternated 

8 

as  in  the  second  two  layups,  one  finds  that  the  load  carrying  ability  of  the 
laminate  is  significantly  improved.  The  onset  stresses  for  each  of  the  three 
fracture  modes  in  the  latter  two  layups  (c)  and  (d)  in  Fig.  10  show  no  signif¬ 
icant  differences,  yet  the  sequence  (d)  which  provides  a  symmetric  [25/-25/-25/ 
25]  layer  on  each  side  of  the  90  layer  show  significantly  higher  load  carrying 
capacity.  It  is  postulated  that  some  l^ss  of  strength  is  associated  with  the 
twisting  of  the  non-symmetric  layup  [25/-25/25/-25J  in  a  region  where  it  is 
effectively  debonded  from  the  90  degree  ply  because  of  the  extensive  edge  de- 
lamination. 

CONCLUSIONS 

It  is  not  the  purpose  of  this  paper  to  provide  a  detailed  comparison 
between  the  experimental  observations  of  fracture  sequence  and  analytical  models 
for  prediction  of  the  Initiation  and  growth  of  specific  fracture  modes.  The 
references  cited  in  this  paper  provide  Information  on  the  more  successful  ap¬ 
proaches  to  modeling  individual  fracture  modes.  Rather  we  have  concentrated  on 
the  documentation  of  the  interaction  of  fracture  modes  immediately  prior  to 
failure  and  the  differences  in  the  modes  and  their  interaction  within  a  given 
series  o*f  laminates. 


From  the  observed  diversity  of  fracture  mode  type  and  interaction  found 
in  this  series  of  laminates,  it  is  difficult  to  conclude  that  any  analytical 
model  which  accounts  for  fracture  mode  initiation  and  growth  and  the  inter¬ 
actions  of  modes  in  a  given  laminate  can  be  applied  "apriori"  to  predict  the 
tensile  stength  of  a  different  layup  sequence.  The  detailed  modeling  of  fail¬ 
ure  mode  interactions  provides  primarily  a  quantitative  understanding  of  the 
mechanics  of  that  interaction  and  a  confidence  that  a  particular  mode  inter¬ 
action  resulting  in  lower  strength  or  stiffness  will  not  occur. 

With  the  exception  of  the  [25/-25/25/-25/90/901  and  [25/-25/-25/25/90/901 
laminates  and  laminates  containing  very  thick  90  deg.  plies,  the  strength  of 
these  [theta,  -theta,  90]g  laminates  is  found  to  be  very  simply  estimated  by 
ignoring  stress  concentrations  associated  with  transverse  cracking  and  de¬ 
lamination  and  assuming  that,  near  the  UTS,  the  load  will  be  carried  entirely 
by  the  primary  load  carrying  ply.  The  plots  of  ultimate  load  in  Figs.  9  and 
11  bear  this  out. 

If  one  assumes  that  eventually  all  of  the  tensile  load  is  redistributed 
onto  the  primary  load  carrying  plies,  a  lower  bound  on  strength  can  be  estab¬ 
lished.  For  more  complex  laminates  like  the  [0/45/-45/90]g  family  this  assump¬ 
tion  will  be  very  conservative  for  certain  stacking  sequences  in  which  delami¬ 
nation  is  not  observed  under  static  tension  [1].  On  the  other  hand,  this  lower 
bound  estimate  of  strength  was  successfully  utilized  by  Ryder  [12]  to  predict 
the  tens ion- tens ion  fatigue  life  and  strength  of  quasi-is tropic  laminates. 

The  results  described  in  this  paper  reinforce  the  argument  that  some  detailed 
experimental  observations  of  the  particular  fracture  sequence  in  a  given  laminate 
are  clearly  necessary  for  accurate  strength  prediction,  although  the  trends  ob¬ 
served  in  a  similar  laminate  construction  may  be  sufficient  to  place  some  useful 
bounds  on  the  predicted  strength. 
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90-Degree  Ply  Cracking(TC)  Onset  Stress  and  Ultimate  Tensile  Strength  As 
Influenced  by  the  Thickness  of  the  90-Degree  Layer. 


90-Degree  Ply  Cracking(TC)  Onset  Stress  and  Ultimate  Tensile  Strength  As 
Influenced  by  the  Thickness  of  the  0-Degree  Layer. 
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Average  Ultimate  Tensile  Load  As  A  Function  of  Angle  and  Ply  Thickness  In 
( theta/-theta/90  )  T300/934  Laminates. 


ULTIMATE  LOAD  VS  STACKING  SEQUENCE 
In  (2G.-26,00)a  LAMINATES 


T300/934  Laminate  family 
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ON  CRACK  DEVELOPMENT  IN  GRAPHITE  EPOXY 


t°2/9°n] s  LAMINATES  UNDER  UNIAXIAL  TENSION 

A.  S.  D.  Wang,  N.  N.  Kishore  and  C.  A.  Li 
Drexel  University 
Philadelphia,  PA  19104 

ABSTRACT.  When  a  graphite  epoxy  laminate  of  cross  ply  (0/90)  lamination 
are  subjected  to  uniaxial  tension,  transverse  cracks  in  the  90°  ply  can 
be  induced  at  strains  as  low  as  0.3Z.  As  the  applied  tension  increases, 
the  cracks  continue  to  form  in  increasing  numbers. 

At  some  late  stage  of  loading  and  before  the  breaking  of  fibers  in 
the  0®  ply,  several  other  matrix  dominant  cracking  modes  may  emerge,  all 
stemming  from  the  sites  of  .the  transverse  cracks.  Among  them  are 
(1)  splitting  cracks  in  the  0s  ply  that  cross  over  one  or  more  90s  ply 
transverse  cracks;  (2)  localized  delamination  in  the  0°/90®  interface 
emanating  from  the  point  where  a  0*  ply  split  crosses  a  90°  ply  transverse 
crack;  and  (3)  localized  edge  delamination  in  the  0®/90®  interface  emana¬ 
ting  from  the  corner  point  where  a  transverse  crack  meets  the  free  edge  of 
the  laminate  coupon. 

These  late  stage  matrix  cracking  modes  have  been  observed  to  form 
under  either  static  or  cyclic  loads.  Their  individual  growth  and  mutual 
coalescence  often  precipitate  disintegration  of  the  lamination  structure 
and  eventuate  its  total  failure. 

The  physical  mechanisms  and  the  load-dependent  crack  development 
process  of  the  90*  ply  transverse  cracks  have  previously  been  described  by 
a  stochastic  model  that  combines  the  general  concept  of  fracture  mechanics 
and  a  random  distribution  of  material  flaws. 


The  objective  of  this  paper  is  to  apply  the  previously  developed 
approach  and  describe  the  various  cracking  events  in  cross-ply  laminates 
as  they  occur  subsequent  to  the  formation  of  90°  ply  transverse  cracks. 

The  particular  laminates  modelled  are  the  [^2/90q]8>  n  ■  1,  2  and  4 
family,  subjected  to  static  tensile  loading.  Emphasis  is  placed  on  de¬ 
termining  why  a  certain  mode  of  cracking  occurs  while  other  modes  do  not. 
Especially  for  the  cracking  state  which  exists  during  late  stage  loading, 
an  analysis  will  be  made  to  explain  why  that  state  is  critical  to  the  final 
strength  of  the  laminate. 

All  predicted  results  will  be  compared  with  the  experimentally 
obtained  data.  The  details  of  the  experiment  have  been  reported  in  an 
earlier  paper. 


I.  INTRODUCTION 


Structural  laminates  made  of  unidirectional  fiber  reinforced  plies 
frequently  suffer  interply  and  intraply  matrix  cranks  well  before  the  final 
rupture  of  the  laminates  QlU .  Consider  as  an  example  the  graphite-epoxy 
cross-ply  (0/90)  type  laminates  subjected  to  uniaxial  tension.  Matrix 
cracks  in  the  90°  ply  known  as  transverse  cracks  have  been  observed  to 
occur  at  laminate  tensile  strains  as  low  as  0.2  to  0.3£  iz, 3] .  Since  the 
final  rupture  of  the  laminate  is  primarily  controlled  by  the  0°  plies,  whose 
ultimate  tensile  strain  exceeds  1 %,  matrix  cracking  of  this  kind  has  often 
been  tolerated  in  design  as  some  benign  internal  failure  mode. 

Actually,  it  has  been  known  that  the  formation  of  intraply  matrix  cracks 
signifies  the  beginning  of  a  complex  sublaminate  crack  development  process, 
which  almost  always  leads  to  cracks  in  multiple  plies  and  internal  load  re¬ 
distribution.  In  the  case  of  the  cross-ply  laminate  cited  above,  onset  of 
90°  ply  transverse  cracking  is  followed  by  multiple  cracks  in  that  ply  if 
the  applied  tension  is  increased.  Or,  if  it  is  cyclically  applied,  multiple 
cracks  in  that  ply  will  form  in  time.  In  either  case,  the  presence  of  a 
distribution  of  the  cracks  in  the  laminate  can  reduce  the  threshold  condition 
for  other  matrix  crack  modes  to  emerge  D*»5.6,?3*  One  such  mode,  for  instance 
is  the  intraply  splitting  of  the  0°  ply,  the  split  being  parallel  to  the 
applied  tension.  Under  ascending  load,  or  in  time,  0°  ply  splitting  also 
occur  in  multiple  numbers.  At  each  crossing  point  where  a  0°  ply  split  inter¬ 
sects  a  90°  ply  transverse  crack,  interaction  between  the  two  crossing  cracks 
causes  a  localized  but  highly  magnified  three  dimensional  stress  field.  The 
tensile  and  the  shearing  components  in  this  stress  field  may  in  turn  cause 
0/90  interface  delamination  around  the  crossing  point  Coalescence  of 

the  many  localized  delaminations  eventuates  a  massive  0/90  interface  separa¬ 
tion  at  some  critical  loading  level.  This  is  often  the  last  matrix  cracking 
event  before  the  catastrophic  failure  of  the  0°  ply  by  fiber  breaking  £5], 


Laminates  of  more  practical  lay-up  suffer  similar  matrix  crack  develop¬ 
ment,  though  the  sequence  of  events  may  not  necessarily  follow  the  same  order, 
Experiments  have  generally  established  that  the  entire  process,  including 
the  cracking  state  which  exists  just  prior  to  final  laminate  failure,  is 
a  generic  phenomenon  unique  to  the  laminate  type  and  the  nature  of  the  load; 
the  phenomenon  may  be  amenable  to  some  form  of  analytical  description  £4]. 

While  such  a  prospect  remains  largely  unreached  at  the  present  time, 
several  successful  efforts  have  been  made  to  describe  crack  development  in 
laminates  of  simpler  forms.  The  basic  mechanisms  of  90°  ply  transverse  crack 
and  free  edge  induced  delamination  have  been  explained  and  predicted  by  the 
energy  release  rate  concept  of  the  classical  fracture  mechanics  £8,9.10]. 
Multiple  cracks  in  the  90°  ply  are  described  by  a  stochastic  procedure  in¬ 
corporating  a  new  concept  of  effective  flaw  distribution  as  an  inherent  ply 
material  property  £ll].  The  utility  of  the  same  approach  to  predict  the 
growth  behavior  of  free  edge  delamination  under  fatigue  loading  has  also 
been  established  £l2]. 

However,  matrix  cracks  in  multiple  plies  and  their  mutual  interactions 
have  scarcely  been  analytically  approached.  In  particular,  an  undersatandlng 
is  needed  regarding  the  crack  state  which  exists  just  prior  to  final  failue 
of  the  laminate.  (Questions  pertaining  to  why  and  how  that  state  has  come 
about,  and  why  is  it  critical  to  the  laminate's  final  failure  axe  yet  to  be 
answered. 

In  an  effort  to  understand  the  mechanisms  of  matrix  crack  interactions, 
a  research  has  been  conducted  which  included  an  experimental  investigation 
and  an  analytical  simulation. 

The  experiment  detailed  the  chronological  matrix  cracking  events  in  a 
family  of  laminate  coupons  in  the  form  of  £±9/90  ].  The  coupons  were  made 
of  the  tJ00/93k  graphite-epoxy  composite  system,  and  were  tested  under  static 


uniaxial  tension.  The  outer  angle-ply  orientation  ©  varied  from  0  ,  25  ,  30 
and  60°{  the  number  of  the  90°  plies  in  the  core  varied  from  n  =  1,2,4  and  8. 
The  ranges  of  these  lamination  variables  were  so  chosen  because  they  effect 
different  matrix  cracking  events  at  distinctly  different  threshold  stress 
conditions.  All  cracking  patterns  were  detected  and  documented  throughout 
the  entire  course  of  loading. 

Results  from  this  experimental  investigation  have  been  presented  in  an 
accompanying  paper  £13!)  • 

This  paper  presents  the  results  from  the  analatical  investigation.  The 
energy  release  rate  method  developed  earlier  for  modelling  the  90°  ply 
transverse  cracks  and  free  edge  delaminations  was  employed  again  to  model 
matrix  crack  interactions.  The  stress  fields  near  several  important  crack- 
interaction  regions  were  analyzed  first  by  a  three  dimensional  finite  element 
routine,  and  the  same  routine  was  used  to  simulate  any  possible  crack  growth 
stemming  from  the  interacting  cracks.  Special  attention  was  placed  on  the 
various  material  and  geometrical  factors  which  play  dominant  role  in  the 
interaction  effects. 

For  conciseness,  however,  only  the  results  for  Co_/90  1  ,  n  «  1,2  and 
4  will  be  reported  here.  As  it  may  become  clear  later,  the  same  method  of 
analysis  can  be  equally  applied  to  laminates  of  any  lay-up  form. 


II.  CRACK  DEVELOPMENT  IN  [Cu/90  1 

c,  ns 

Let  us  consider  the  laminate  family  £oo/90  1  ,  for  n  =  1,2  and  4.  Coupons 

n  s 

of  25  mm  (1  in.)  width  and  225  mm  (9  in.)  length  were  made  and  tested  under 
monotonic  tension  in  room  temperature  condition.  The  material  properties  of 
the  T 300/934  unidirectional  ply  have  been  previously  characterized  and  are 
listed  here  in  Table  I.  As  has  been  mentioned,  the  details  of  the  test  and 
test  results  are  reported  in  Ref.  £l33*  and  the  energy  method  used  for  model¬ 
ling  is  found  in  Refs.  C8,9»103«  The  three  dimensional  finite  element  routine, 
it's  development  and  user’s  instructions  are  contained  in  Ref.  £l43.  In  this 
section,  we  simply  apply  these  techniques  and  predict  the  sequential  crack 
events  as  the  laminates  are  loaded  in  monotonic  tension. 

Transverse  Cracking  Versus  Free  Edge  Delamination. 

When  the  COg/^^g  laminates  are  loaded,  two  competing  crack  modes  exist 
during  the  initial  stage  of  loading.  One  mode  is  due  to  the  in-plane  stress  In 
the  90°  ply  which  causes  multiple  cracks  in  that  ply  as  shown  schematically 
in  Fig.  1.  The  other  possible  cracking  mode  is  free  edge  delamination  caused 
by  interlaminar  tensile  stress  which  is  concentrated  near  the  free  edge  £15, 
163.  A  schematics!  view  of  the  latter  is  shown  in  Pig.  2. 

It  will  be  determined  next  that  the  free  edge  effect  in  the  £02/90n3g 
family  under  tension  is  negligible;  and  edge  delamination  alone  cannot  occur. 
Consequently,  90°  ply  transverse  cracking  is  the  first  failure  event. 

To  show  this,  we  first  calculate  the  energy  release  rate  coefficient 
functions  Cg(a)  and  C,p(a)  associated  with  90°  ply  transverse  cracking.  The 
significance  of  these  coefficients  is  that  the  available  energy  release  rate 
G  for  a  flaw  in  the  90°  ply  to  become  a  transverse  crack  can  be  expressed  as 


G(i  ,  AT, a)  =  [>/£  -5V  +  VV  AT  ]2  t 


(1) 


where  ex  is  the  applied  laminate  tensile  strain,  AT  is  the  temperature  load 
experienced  by  the  laminate  after  curing  and  t  is  a  length  parameter  taken 
here  as  one  ply  thickness.  The  quantity  a  represents  the  size  of  the  assumed 
effective  flaws  (=2a)  which  distribute  in  the  90°  plies(  for  details  see  ClcQ 
or  [ll]) . 

Figures  3  and  4  show  the  calculated  Cg( a)  and  C^(a)  functions  for  the 
laminates  and  COg/^O^g*  Using  these  curves,  the  threshold  stress 

condition  for  the  first  transverse  crack  and  the  subsequent  cracks  can  be 
obtained.  This  will  be  shown  later  in  this  section. 

If  we  assume  that  free  edge  delamination  will  occur  before  transverse 
cracking,  we  then  calculate  the  energy  release  rate  coefficient  functions  C 

e 

and  associated  with  an  edge  flaw  of  size  a(for  details  see  Qlol  and  £12^) . 
Figures  5  and  6  show  the  calculated  C  (a)  and  C_(a)  curves  for  Qoo/90„1  and 
Co2/904]s  assuming  the  edge  flaw  in  the  laminate  mid-plane  causes  delamination. 
These  curves  can  be  used  to  predict  the  onset  load  condition  for  delamination. 

The  criterion  for  the  flaw  to  become  a  crack  is 

G(Sx,AT,a)  =  Gc  (2) 

where  Gq  is  the  appropriately  determined  critical  energy  release  rate  of  the 
materlal(see  discussions  in  £173) • 

Onset  of  Transverse  Cracking.  To  determine  the  onset  load  for  90°  ply 
transverse  cracking,  the  size  of  the  most  dominant  effective  flaw,  2aQ,  in  the 
effective  flaw  distribution  must  be  identified.  Following  the  reasoning  given 
in  Ref.  is  defined  as  the  99  percentile  of  the  effective  flaw  dis¬ 

tribution  f(a).  As  it  will  be  shown  later,  a  for  both  Lo«/907]  and  Co,/90k3 
is  in  the  order  of  1.5t.  Thus,  from  Figs.  3  and  4  we  obtain, 
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for  n  =  2, 


for  n  =  4, 


C  (a  )  *  41.3x10''  J/nr 
e  0  A 

(6x10°  in-lb/in-5) 

CT(ao)  =26.1  j/m3/C2 

(11. 7x10 ~4  in-lb/in^/F2) 

Ce(aQ)  =  48.3xl09  J/m3 

(7xl06  in-lb/in3) 

CT(ao)  =26.05  J/m3/C2 

(ll.7xl0“4  in-lb/in3/?2) 


Substitution  of  these  values  into  Eq.,(l)  and  using  AT  =  125°C  (225°F), 

t  =  0.132  mm  (0.0052  in.)  and  G  =  228  J/m2  (1.3  in-lb/ln2)  from  Table  I,  we 

c 

obtain  from  Eq..(2)  the  respective  onset  strains  for  transverse  crackings 


for  n  =  2,  (Sx)tc  =  3-3  x  10~3 

for  n  =  4,  (ex)tc  =  3.1  x  10~3. 


Onset  of  Free  Edge  Delaml nation.  If  we  assume  free  edge  delamination  in 
the  mid-plane  to  occur  before  transverse  cracking,  the  associated  threshold 
stress  can  be  determined  by  using  the  maximum  values  of  the  Cg  and  curves 
in  Figs.  5  and  6.  In  essence,  we  assume  the  size  of  the  effective  edge  flaw 
aQ  >  2t;  and  the  prediction  yields  the  lower  bound.  Hence, 


for  n  =  2, 


C 

e,max 

CT,max 


=  0.73  x  109  J/m3 
(0.106  x  106  in-lb/in3) 

=7-03  J/m3/C2 
(3.15  x  10-4  in-lb/in3/?2) 


for  n  =  4,  C  =  0.41  x  109  J/m3 

e.max  ,  '  - 

(0.06  x  10  in-lb/in-5) 

°T,max3  V»3/C2 

(1.67  x  10-4  in-lb/in^F2) 


Using  the  sane  values  for  AT  and  t  as  before,  but  G  =  158  J/ni  (0.9 

c 

2 

in-lb/in  )  for  mode-I  mid-plane  delamination  as  indicated  in  Table  I,  we 
obtain  from  Bq,s.  (1)  and  (2)  the  respective  onset  strains( lower  bounds)  for 
free  edge  delamination : 

for  n  =  2,  (®x)de  =  36  x  10"3 

for  n  *«  4,  (®x)de  =  52  x  10~3* 

Clearly,  the  assumed  free  edge  delamination  will  never  occur,  because 
the  laminate  final  failure  by  fiber  breaking  is  likely  to  take  place  at 
about  10  x  lo'3  to  12  x  10-3  strain. 

Influence  of  9  in  £±9/90]  Laminates.  The  previous  examples  showed  that 

s 

free  edge  effect  is  relatively  unimportant.  However,  when  the  angle  9  increases 

slightly  from  being  0°,  the  role  of  free  edge  effect  increases  dramatically. 

Fig.  7  shows  a  family  of  the  C  (a)  curves  associated  with  delamination  in 

© 

the  laminate  mid-plane  for  C±9/90]  with  9  =  0°,10°t20o  and  30°.  It  is  seen 

s 

that  Cg  rises  sharply  with  9.  On  the  other  hand,  Fig.  8  shows  the  correspond¬ 
ing  family  of  the  C^,  curves,  which  decrease  with  9.  Between  these  two  sets  of 
curves,  the  angle  9  =  25°  yields  approximately  the  optimum  free  edge  effect. 

In  fact,  if  we  use  the  same  calculation  procedure  as  before,  the  onset  strains 
for  mid-plane  delamination  in  Co,/90]  and  in  C ±25/90]  are: 

for[02/90]s,  (5x)de  =  24  x  10_3 

for[±25/90]a,  ( ®x) de  =  6<  x  10’3* 

From  these  results,  it  may  be  concluded  that  free  edge  delamination  is  not 
possible  for  the  case  of  9  =  0°,  but  is  definitely  possible  for  the  case  of 
9  =  25°.  It  is  not  known  however,  whether  transverse  cracks  in  the  90°  ply 


occur  before  or  after  free  edge  delamination.  So,  let  us  Investigate  the 

Influence  of  9  on  transverse  cracking.  Fig.  9  shows  the  C  (a)  family  associated 

0 

with  a  transverse  crack;  and  the  corresponding  C,p(a)  family  is  displayed  in 

Fig.  10.  It  is  seen  that  6,  in  the  range  indicated,  has  only  a  small  effect 

on  transverse  cracking.  In  fact,  the  calculated  onset  strains  for  Co2/90]g 

and  £ ±25/90l  are i 
s 

for  [02/90]s,  (ix)tc  =  5* 9  x  10 "3 

for  [±25/90]s,  (5x)te  =  6,3  x  10~3‘ 

From  the  above  results,  we  can  readily  determine  that  transverse  cracking 

in  C+25/90!  laminate  should  occur  after  free  edge  delamination, 
s 

Summary.  In  this  section,  we  have  demonstrated  the  profound  influence  of 
lamination  variables  on  the  formation  mechanisms  of  two  competing  matrix  crack 
modes;  namely,  transverse  cracking  and  delamination.  The  energy  release  rate 
method  takes  into  account  these  variables  and  provides  a  quantitative  predic¬ 
tion.  Table  II  summarizes  the  predicted  results  obtained  in  this  section;  and 
these  are  compared  with  the  corresponding  experimental  results  extracted  from 
Ref.  [13]. 

90°  Ply  Multiple  Cracks  in  Co9/90  !  . 

Ct  ns 

In  the  previous  section,  it  was  shown  that  free  edge  delamination  as  a 
first  matrix  failure  mode  will  never  occur  in  the  C0-/90  !  family  when  they 
are  loaded  in  tension.  Rather,  the  first  matrix  failure  mode  is  the  90°  ply 
transverse  cracking  which  Initiates  at  distinctively  different  laminate  tensile 
strains,  depending  on  the  value  of  n  or  the  number  of  the  90°  plies.  Following 
the  onset  of  the  first  crack  in  the  90°  plies,  multiple  cracks  can  occur  as 


the  applied  tension  Increases,  resulting  in  a  characteristic  crack  density 
versus  applied  load  relationship.  The  stochatlc  search  procedure  based  on  an 
assumed  effective  flaw  distribution  can  be  applied  to  simulate  the  crack  de- 


velopment  process.  The  procedure  was  previously  formulated  in  Ref.  £lOj  and 
was  discussed  In  more  details  In  Ref.  LllD.  We  shall  therefore  omit  the  compu 
tatlonal  details  here  and  simply  present  the  simulated  results. 

The  simulation  will  be  based  on  the  assumed  flaw  distribution  which  is 
characterized  by  the  flaw  size  distribution  and  the  flaw  spacing  distribution 
These  two  distributions  are  represented  by  the  probability  density  functions. 


f(a)  =  -Jr" 
a>/55f 

exp  [-(a  -  4a)2/2Va2] 

(3) 

1 

f(s)  = 

exp  [-(s  -  4  )2/2v  2] 

(4) 

aV2rr 

8  S 

where  4  and  v  denote  the  mean  and  the  standard  deviation  of  the  normal  dis- 
trlvution  function  f(). 

For  the  T300/934  material  system  used  here,  the  parameters  4  and  v  in 
the  assumed  flaw  distribution  functions  were  previously  determined  as  (see 
Ref.  DU)]). 

for  n  =  1,  4  =  0.053  mm  v  =  0.019  mm 

0.0021  in.  0.00075  in. 

4  =  0.313  mm  V  =  0.115  mm 

s  s 

0.0125  in.  0.00*46  in. 

for  n  =  2,4,  4=  0.091  mm  v  =  0.033  mm 

&  & 

0.0036  in.  0.0013  in. 

4  =  0.313  mm  v  =  0.115  mm 

s  s 

0.0125  in.  0.00*16  in. 

Now,  let  the  size  of  the  roost  dominant  flaw  In  each  case  be  denoted  by 
2aQ  and  let  aQ  be  the  99  percentile  of  f(a).  Then,  for  the  normal  distribu¬ 
tion  expressed  In  (3),  aQ  can  be  approximated  by  £183, 

%  =  +  3Va  (5) 


It  follows  that  the  largest  flaw  determines  the  onset  of  transverse 


cracking,  and  their  size  is  found  using  (5).  Hence,  for  n  =  1,  aQs  0.112  mm 
(0.0044  in.) j  for  n  =  2  and  4,  ao=  0.191  mm  (0.0075  in.).  The  former  is  about 
0.8  times  the  ply  thickness  t,  while  the  latter  is  about  1.5t. 

The  stochastic  procedure  actually  mimics  the  cracking  process  as  it 
would  occur  naturally.  The  result  is  a  simulated  crack  density  versus  applied 
load  relationship.  Figures  11,12  and  13  display  such  a  relationship  for,  res¬ 
pectively,  the  n  =  1,2  and  4  laminates.  The  simulated  results  were  generated 
for  four  or  five  specimens  in  each  case;  they  collectively  form  a  scatter 
bend  as  shown  in  each  of  these  figures.  The  corresponding  experimental  result 
is  shown  by  the  shaded  bend. 

Of  course,  the  simulation  assumes  that  the  applied  tension  can  be  mono- 
tonically  increasing;  and  the  resulting  transverse  cracking  process  can  go 
on  forever.  This,  clearly,  does  not  happen  in  reality  since  other  modes  of 
matrix  cracking  may  emerge  at  higher  load  levels.  As  it  will  be  shown  next, 
any  late-stage  matrix  cracking  invariably  involves  interactions  with  the 
transverse  cracks  in  the  90°  plies.  These  distributed  cracks  act  like  local 
defects,  each  is  capable  of  causing  localized  stress  concentration. 

Matrix  Cracks  At  High  Load  Levels. 

When  a  transverse  crack  is  formed  in  the  90°  ply,  the  load  formerly 
carried  by  the  90°  ply  must  now  be  taken  by  the  0°  ply  near  the  vicinity  of 
the  transverse  crack.  This  extra  load,  however,  is  gradually  transferred  back 
to  the  unbroken  part  of  the  90°  ply  through  the  bonding  strength  of  the  0/90 
interface.  Thus,  near  the  transverse  crack  root  ,  a  stress  concentration 
zone  exists  which  may  cause  both  in-plane  and  out-of-plane  matrix  cracking 
modes . 


In  addition,  the  presence  of  free  edges  now  becomes  Important  factor. 

As  is  shown  schematically  in  Figure  14,  a  highly  concentrated  stress  zone 
exists  near  the  intersection  of  a  transverse  crack  and  the  coupon  free  edge. 


We  shall  label  this  point  as  the  'corner  point'.  The  high  stresses  near  this 
point  are  caused  by  the  compound  effects  from  the  transverse  crack  and  the 
free  edge. 

The  effects  of  these  localized  stress  concentrations  on  matrix  cracking 
at  high  loads  will  be  investigated  in  this  section.  Since  these  local  stress 
fields  are  three  dimensional,  the  finite  element  routine  developed  previously 
Cl4]  will  be  employed. 

Stress  Field  Wear  A  Transverse  Crack.  The  tensile  coupon  shown  in  Fig. 14 
already  has  multiple  transverse  cracks  in  the  90°  ply.  For  purpose  of  analysis, 
consider  only  a  section  of  the  laminate  which  contains  one  crack.  Due  to 
lamination  symmetry,  it  Is  sufficient  to  model  one-eighth  of  this  section  by 
a  net-work  of  finite  elements.  The  actual  sizing,  meshing  and  other  related 
computational  details  are  found  in  Ref.  £l4].  Here,  we  shall  discuss  primarily 
the  stresses  which  act  on  the  0/90  interface  near  the  crack  root,  and  the 
stresses  acting  in  the  0°  ply  above  the  crack  root.  It  is  noted  that  the  inter¬ 
face  stresses  can  cause  delamination,  while  the  in-plane  stresses  can  cause 
several  possible  failure  modes  in  the  0°  ply. 

Figures  15,  16  and  17  show,  respectively,  the  Interlaminar  stresses  d^, 

Txz  and  T  which  are  acting  on  the  0/90  interface.  Figures  18,  19  and  20  show 

the  in-plane  stresses  d  ,  d  and  t  ,  respectively,  which  are  acting  in  the  0° 

x  y  xy 

ply  but  close  to  the  0/90  interface.  Both  sets  of  figures  are  for  the  cases 
of  n  =  2  and  4.  And,  the  magnitude  of  the  stresses  is  due  to  a  unit  tensile 
strain  applied  at  the  far- field,  ex=  10 

Similar  display  for  the  interlaminar  and  the  in-plane  stresses  due  to  a 
unit  of  temperature  change  is  shown  in  Figures  21  to  26,  the  unit  temperature 
being  AT  =-0.56°C  (-1°F) . 

In  all  the  stress  displays,  the  xy -plane  coincides  with  the  0/90  interface, 


with  the  x-ajds  lying  along  the  free  edge  of  the  laminate,  and  the  y-axis 
lying  along  the  transverse  crack  root)  the  origin  of  the  axes  is  the  corner 
point  shown  in  Fig.  14. 

A  cursory  examination  of  these  stress  displays  reveals  the  foilwing 
general  features! 

1.  The  influence  of  laminate  free  edge  on  the  interlaminar  and  the  in-plane 
stresses  is  relatively  insignificant  compared  to  the  influence  of  the  90° 
ply  transverse  crack.  But  the  two  effects  compound  each  other  neax  the 
comer  point,  resulting  in  a  much  amplified  stress  concentration  zone. 

The  amplification  of  stress  near  this  point  is  not  a  simple  addition  of 
these  two  individual  effects. 

2.  The  influence  of  the  90°  ply  thickness  (the  value  of  n)  on  the  stresses 
continues  to  dominate.  An  increase  of  n  increases  both  the  magnitudes  of 
the  stresses  and  the  size  of  the  stress  concentration  zone.  Clearly,  any 
material  damage  Inside  this  zone  will  be  likewise  controlled  by  the  same 
factor. 

3.  The  stress  fields  caused,  respectively,  by  ex=  lO-^  and  by  AT  =  -0.56°C. 
have  their  corresponding  stress  components  in  the  same  sign  inside  the 
concentration  zone.  For  instance,  due  to  tension  is  positive  neax  the 
corner  point  and  along  the  transverse  crack  root  (Fig.  15) t  and  O'  due  to 
the  temperature  load  is  also  positive  in  the  same  region  (Fig.  21).  Thus, 
they  combine  to  give  a  stress  field  of  much  higher  magnitude  in  the  actual 
laminate. 

It  will  be  shown  in  the  sequel  that  the  local  stress  concentration  which 
surrounds  a  transverse  crack  in  the  90°  ply  is  capable  of  causing  different 
modes  of  material  damage  at  high  loads.  And  the  thickness  factor  (the  value 
of  n)  plays  an  important  role  in  this  regatrd . 


•„v 


Effects  of  the  Interlaminar  Stresses  -  Localized  Delamlnatlon.  As  has 


been  shown  in  Figures  15  and  21.  the  normal  stress  (?z  acting  on  the  0/90  inter¬ 


face  is  tensile  near  the  corner  point  and  along  the  transverse  crack  root. 

This  stress  alone  can  debond  the  interface  at  least  locally.  In  addition,  the 
shearing  stresses  r ^  and  t  are  also  of  considerable  magnitude  (see  Figures 
16,  17  and  22,  23).  Hence,  the  actual  debonding  will  be  essentially  mixed-mode 
in  nature.  Vhen  the  applied  load  reaches  a  certain  magnitude,  a  contour-shape 
delamination  will  occur  emanating  from  the  comer  point.  This  is  schematically 
illustrated  in  Figure  27. 

In  order  to  make  a  prediction  for  the  onset  and  the  growth  behavior  of 
the  contour  delamination,  we  begin  with  the  assumption  that  a  small  interface 
flaw  is  located  near  the  comer  point.  This  assumption  is,  of  course,  an  out¬ 
growth  of  the  effective  flaw  concept  which  was  discussed  earlier  in  this  paper. 
Under  the  applied  load,  the  assumed  flaw  will  become  a  local  delamination  crack 
having  a  contour  shape.  The  growth  of  the  contour  is  mimicked  by  a  sequence 
of  nodal  point  release  in  the  finite  element  net-work.  This  is  represented  by 
an  approximate  sequential  nodal  release  as  displayed  in  Figure  27.  As  each  of 
the  finite  element  nodes  is  released,  the  corresponding  .delaminated  area  as 
well  as  the  strain  energy  lost  can  be  computed  (for  details,  see  Ref.  Cl4]) . 

In  this  manner,  a  relationship  between  the  delaminated  area  and  the  strain 
energy  release  rate  is  obtained  for  the  assumed  contour-shaped  growth. 

Figure  28  shows  the  computed  energy  release  rate  coefficient  Cg( for  1) 
as  a  function  of  the  assumed  delamination  contour  for  the  Co,/90j  laminate. 
Note  that  Cg  consists  of  all  three  modal  components,  each  of  which  is  also 
shown  in  the  figure.  In  this  case,  the  mode-I  component  is  slightly  larger 
than  the  mode-II  component,  while  the  mode-III  component  is  negligible.  Also, 
the  release  of  those  nodes  along  the  transverse  crack  root  yields  the  largest 
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energy  release  rate. 

Figure  29  shows  the  coefficient  for  the  laminate  subjected  to  AT  =  -1*F 
or  -0.65*C.  The  general  character  of  C^,  and  it's  three  modal  components  are 
similar  to  those  of  C  . 

Although  the  calculated  C  and  C_  coefficients  axe  based  on  the  assumed 

©  A 

nodal  release  sequence  (shown  in  Figure  27) t  which  may  not  represent  the  best 

possible  choice,  the  general  trend  of  these  coefficient  functions  indicates 

that  the  delamination  growth  is  basically  a  stable  process.  This  qualitative 

result  agrees  with  the  actual  behavior  observed  in  experiment. 

Furthermore,  if  we  assume  that  the  actual  delamination  becomes  macro- 

2 

scopically  observable  when  it's  area  grows  to  about  t  (t  is  one  ply  thickness), 

we  then  define  such  a  state  as  the  onset  of  the  contour  del anl nation.  Thus,  we 

take  the  maximum  values  of  C@  and  C^,  from  Figures  28  and  29  at  k/\T  =  1  and 

compute  the  total  energy  release  rate  G(ex,AT,A)  via  Eq.  (l) .  Let,  as  before, 

AT  =  -125°C  (-225°F)  and  G  =  289  J/m2  (1.65  in-lb/in2)  for  mixed-mode  crack 

c 

growth,  see  Table  I.  It  follows  from  Eq.  (2)  that  the  threshold  tensile  strain 
for  the  onset  of  the  contour  delamination  is  predicted  as 

(ix)de  *  11.0  x  10 "3  (for  n  *  2)  (6) 

Similarly,  the  computed  Cg  and  coefficient  functions  for  are 

shown  in  Figures  30  and  31 »  respectively.  It  is  seen  that  the  values  of  Gg  and 
CT  are  significantly  increased  as  a  result  of  doubling  the  number  of  the  90° 
plies.  Consequently,  the  predicted  threshold  strain  for  onset  of  the  contour 
delamination  is  much  reduced | 

(•x>de  =  7,8  x  10"3  (for  n  =  4)  (?) 

Now,  if  we  consider  the  fact  that  the  breaking  strain  of  \he  0°  ply  is 


10  to  12  x  10  J ,  it  becomes  clear  that  the  assumed  contour  delamination  near 
the  transverse  crack  will  occur  in  the  COg/^^]  laminate  well  before  any  fiber 
break  in  the  0°  ply,  though  it  is  not  likely  to  occur  at  all  in  the  jj^/^JOgDg 
laminate. 

It  is  also  clear  fiom  the  results  in  (6)  and  (7)  that  such  a  type  of 
delamination  will  never  occur  in  the  £o„/90]  laminate  before  the  final  laminate 

c  8 

failure. 


Effects  of  the  In-plane  Stesses  -  Fiber  Break  and  0  Ply  Split.  Let  us 

examine  next  possible  failure  modes  effected  by  the  in-plane  stresses  which 

exist  in  the  0°  ply  and  are  magnified  by  the  presence  of  the  transverse  crack. 

These  in-plane  stresses  are  shown  in  Figures  18,  19  and  20  for  tensile  load  of 

ex=  10"6,  and  in  Figures  24,  25  and  26  for  thermal  load  AT  =-0.65°C  (-1°F). 

Note  that  in  both  cases,  the  normal  stresses,  d  and  d  ,  axe  tensile  and  axe 

x  y 

much  magnified  near  the  transverse  crack;  but  the  shear  stress,  T  ,  is  of 

xy 

small  amplitude  except  near  the  corner  point.  It  will  be  shown  here  that  d^ 
may  cause  local  fiber  break  while  d  causes  0°  ply  splitting  cracks. 

y 

First,  consider  d which  is  taken  by  the  fibers.  The  part  due  to  the 
applied  tension  is  shown  in  Figure  18;  its  amplitude  is  seen  to  increase  30% 
to  50%  near  the  transverse  crack  root,  depending  on  the  value  of  n.  The  part 
due  to  the  thermal  load  is  shown  in  Figure  24.  This  stress  is  compressive 
throughout  the  0°  ply  except  near  the  transverse  crack  root  where  it  becomes 
tensile.  Again,  the  magnitude  of  the  localized  tensile  stress  depends  also  on 
the  value  of  n.  Hence,  in  the  vicinity  of  the  transverse  crack,  the  combined 
<*x  in  the  0°  ply  could  be  as  much  as  50%  to  100%  higher  than  the  far-field 
value.  As  illustrated  in  Figure  32,  d exhibits  a  singular  behavior  approaching 
the  0/90  interface,  and  reduces  magnitude  away  from  the  interface.  It  is  clear 
that  the  localized  concentration  of  the  tensile  stress  will  have  varying  effect 
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on  the  breaking  mechanisms  of  the  0°  ply  fibers,  depending  on  the  90°  ply 
thickness.  But  the  exact  nature  of  fiber  breaking  mechanisms  in  laminates  is 
not  the  object  of  this  paper,  and  we  shall  not  attempt  to  make  any  prediction 
based  on  the  calculated  d*  here.  It  is  sufficient  to  note  that  experimental 
observations  £6,133  have  verified  fiber  breakage  near  the  site  of  90°  ply 
transverse  cracks. 

Consider  next  d  which  is  normal  to  the  fibers  of  the  0°  ply,  see  Figures 

y 

19  and  25.  In  this  case,  both  the  mechanical  part  and  the  thermal  part  axe 
tensile  throughout  the  0°  ply,  but  are  especially  amplified  near  the  transverse 


face  as  is  illustrated  in  Figure  33«  And,  depending  on  the  value  of  n,  the 
maximum  amplification  can  be  as  much  as  more  than  500%.  Such  a  high  stress  is 
capable  of  causing  matrix  cracks  in  the  interior  of  the  0°  ply,  known  as 
longitudinal  splitting. 

The  actual  mechanisms  of  0°  ply  splitting  are  similar  to  that  in  the  90° 
ply  transverse  cracking.  Hence,  the  crack-closure  precedure  and  the  criterion 
used  for  defining  the  onset  of  transverse  cracking  can  be  applied  to  0°  ply 
splitting.  But  we  shall  omit  the  details  of  the  calculation,  and  present  here 
simply  the  simulated  results.  Namely,  matrix  splits  in  the  interior  of  the  0° 
ply  will  occur  over  the  transverse  crack  when  the  applied  tensile  strain  ex 
reaches  the  following  critical  value: 


for  n  =  2, 

(5x)sp  -  10.5  x  io’3 

(8) 

for  n,=  4, 

*  7-8  x  10'3 

(9) 

According  to  the  above  prediction,  0°  ply  splitting  may  or  may  not  occur 
1°  when  under  static  tensionj  but  will  definitely  occur  in  Co2/9o4]s 


before  fiber  breaking.  Note  that  in  the  case  of  the  threshold 

condition  for  0°  ply  splitting  is  about  the  same  as  the  threshold  condition 
for  transverse  crack  induced  0/90  Interface  delamination  (see  Eq.7).  Thus, 
these  two  types  of  matrix  cracks  may  actually  occur  at  the  same  load  level 
for  this  particular  laminate. 

Delamination  at  the  Grossing  Cracks.  When  a  splitting  crack  is  formed  in 
the  0°  ply,  it  makes  a  cross  with  the  transverse  crack  in  the  90°  ply,  see 
Figure  34.  In  the  close  vicinity  of  the  crossing  point,  a  severe  stress  con¬ 
centration  is  developed.  Especially  the  interlaminar  stresses  acting  on  the 
0/90  interface  are  capable  of  causing  localized  delamination  as  shown  by  the 
shaded  area  in  Figure  34.  Hence,  we  shall  examine  only  the  distributional 
nature  of  these  interlaminar  stresses. 

Figure  35(a)  shows  the  0/90  interface  plane.  Note  that  the  x-axis  is  made 
parallel  to  the  0°  ply  split  and  the  y-axis  is  parallel  to  the  90°  ply  trans¬ 
verse  crack;  the  origin  of  the  axes  is  the  intersection  point  of  the  two  cross 
cracks.  The  first-quadrant  distribution  of  the  computed  stresses  0^,  T  and 

r  are  displayed,  respectively,  in  Figures  36,  37  and  38.  Again,  these  are 
yz 

for  the  EOg/^OgDg  and  lajnlna'tes  only;  811(1  the  mechanical  part 

—  —6  o 

of  the  stresses  is  due  to  e^  =  10  while  the  thermal  part  is  due  to  AT  =  -1  F 

or  -0.56°C.  It  is  seen  that  all  the  stresses  are  highly  amplified  near  the 

crossing  point;  and  is  tensile  (both  the  mechanical  and  the  thermal  parts) 

z 

around  this  point. 

In  order  to  determine  the  threshold  condition  for  delamination  near  the 


crossing  point,  a  crack  growth  simulation  is  performed  using  the  finite  ele¬ 
ment  routine.  The  simulated  delamination  followed  the  nodal  release  sequence 
shown  in  Figure  35(8).  The  point-by-point  energy  release  rate  coefficients  Cg 
and  C^,  are  shown  in  Figures  39  and  40  for  and  in  Figures  41  and  42 


From  the  energy  release  rate  coefficient  curves,  it  is  apparent  that  the 
behavior  of  the  delamination  is  in  mixed-mode,  mode-I  and  mode-II  being  about 
eq.ual  and  mode-III  negligible.  The  total  energy  release  rate  is  slightly  de¬ 
creasing  with  delaminated  area  A.  This  again  indicates  that  the  delamination 
growth  is  essentially  stable. 

If  we  compare  the  four  energy  release  rate  curves  with  those  shown  in 
Figures  28,  29,  30  and  31,  we  see  the  close  similarity  between  this  case  and 
the  one  discussed  earlier  for  delamination  emanating  from  the  transverse 
crack/free  edge  intersection.  However,  the  calculated  energy  release  rate  for 
delamination  at  the  cross  cracks  is  slightly  higher. 

Indeed,  if  we  follow  the  same  procedure  as  before,  we  obtain  a  prediction 
for  the  threshold  condition.  Thus,  onset  of  delamination  occurs  when  the  applied 
strain  reaches  the  following  values: 


for  n  =  2, 

<5x>de  ’  10'9  *  10'3 

(10) 

for  n  =  4, 

( ®x>  de  *  7'5  *  10"3 

(11) 

It  appears  that  delamination  at  the  cross  cracks  occurs  at  about  the 
same  time  when  the  0°  ply  splitting  is  formed.  Physically,  of  course,  the 
splitting  must  occur  before  delaminatlon.  It  is  also  noted  that  for  the 

laminate  the  splitting  and  the  delamination  are  unlikely  to  occur 
before  fiber  break. 

Summary  of  Results.  All  results  predicted  in  this  section  are  summarized 
in  Table  III,  Whenever  possible,  the  predicted  results  are  compared  with  their 
experimental  counterparts.  The  latter  were  actually  extracted  from  Ref.  £13], 

A  cursory  examination  of  Table  III  indicates  readily  that  the  predicted  events 
and  the  associated  threshold  strains  agree  well  with  the  experimentally  obatlned 


results.  It  shold  be  emphasized  that  the  predictions  involved  some  degree  of 
approximation,  both  in  numerical  accuracy  and  in  setting  the  criterion  for  the 
threshold  condition. 

In  spite  of  such  reservations,  it  seems  that  the  energy  release  rate 

concept  continues  to  be  applicable  to  matrix  cracking  simulation  problems  at 

high  load  levels  where  complicated  crack  interactions  take  place.  This  concept, 

however,  could  not  at  least  at  the  present  time  predict  the  failure  of  the 

fibers,  although  the  event  of  fiber  failure  (hence  also  the  final  failure  of 

the  laminate)  is  affected  by  matrix  crack  interactions  at  high  loads.  The 

adversity  of  the  interactions  intensifies  if  the  matrix  cracks  can  become 

large.  For  instance,  among  the  three  laminates  of  £0-/90  3  ,  n  =  1,  2  and  4, 

c  ns 

matrix  cracks  in  n  =  1  are  the  smallest  in  size,  while  in  n  =  4  are  the  largest. 
Thus,  at  least  qualitatively,  we  may  conclude  that  the  threshold  strain  for 
fiber  break  should  be  the  highest  in  n*l  and  lowest  in  n  =  4.  This  observa¬ 
tion  is  reasonably  supported  by  experiment  (see  results  in  the  last  column  of 
Table  III) . 

Finally,  let  us  note  that  the  predicted  threshold  strains  for  (a)  0°  ply 

splitting,  (b)  delamination  at  the  cross  cracks  and  (c)  delamination  at  the 

free  edge/transverse  crack  comer  point  are  practically  all  the  same  for  the 

C0o/90l  laminate  family,  although  for  laminates  with  n  ^  2  these  crack  modes 
ns 

are  not  likely  to  occur  before  fiber  break.  This  observation  has  also  been  ex¬ 
perimentally  supported  in  at  least  one  instance;  Figure  43  shows  a  sequence  of 
x  -  ray  photographs  taken  during  tension  loading  of  a  Co2/9°3^s  laminate.  At 
about  90#  of  the  ultimate  load,  cross  cracks  were  seen  to  form  in  increasing 
numbers  and  delaminations  around  the  cross  cracks  followed  almost  instantly. 

And,  the  coalescence  of  the  numerous  local  delaminations  finally  triggered 
the  failure  of  the  laminate.  An  identical  specimen  loaded  by  the  same  procedure 
developed  instead  delaminations  near  the  free  edge/transverse  crack  corners, 


also  at  about  90#  the  laminate  ultimate  load,  see  Figure  44,  While  the  general 


appearances  of  the  two  crack  events  at  high  load  are  different,  the  actual 
threshold  conditions  separating  them  Me  close.  The  fact  that  such  results  can 
be  predicted  at  all  is  in  itself  surprising. 

In  Table  III,  experimental  results  for  COg/tJOg]  were  listed  for  comparison. 
An  analytical  prediction  using  the  same  method  can  be  readily  performed,  but 
it  was  not  attempted  in  this  paper. 
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III.  CONCLUSIONS 


In  this  paper,  we  have  demonstrated  that  the  development  process  of  matrix 
cracks  and  crack  interactions  at  high  loads  can  be  analytically  described  by  a 
three  dimensional  stress  analysis  based  on  the  method  of  ply  elasticity  and  by 
a  crack  growth  simulation  routine  based  on  the  method  of  fracture  mechanics  in 
conjunction  with  a  new  concept  of  material  "effective"  flaw  distribution.  The 
laminate  family  made  of  T300/934  graphite-epoxy  system  was  used  to 
illustrate  and  to  obtain  specific  results  which  can  be  compared  by  a  parallel 
experiment.  To  conclude  this  paper,  we  wish  to  make  the  following  comments: 

1.  Matrix  cracks  in  off-axis  plies,  the  90°  ply  transverse  cracks  being  examples, 
are  important  precursors  in  the  laminate  damage  accumulation  process.  In 
time,  or  at  high  loads,  these  cracks  eventually  cause  secondary  matrix 
cracks  into  the  load-carrying  plies.  When  this  happens,  it  is  often  the 
beginning  of  the  end  of  the  laminate's  failure  process.  In  this  paper,  we 
have  shown  that  the  matrix  crack  development  process  is  influenced  largely 
by  the  lamination  geometry,  given  the  material  and  the  loading  environment. 
And,  the  geometrical  factors  can  be  explicitly  included  in  the  predictive 
models.  Consequently,  it  appears  that  the  final  laminate  failure  event 
could  also  be  much  influenced  by  the  same  factors.  The  present  and  past 
experiments  [_1,13D  have  demonstrated  convincingly  this  fact.  To  emphasize, 
let  us  examine  briefly  the  experimental  results  Ll33  shown  in  Table  IV, 
where  the  final  failure  (tensile)  strains  of  the  Ci©/90  ]  laminate  family 
(n  =  1,  2,  4  and  8j  9  «  0°,  25°  and  30°)  are  displayed.  It  is  seen  that  the 
failure  strains  are  influenced  by  the  value  of  n.  And,  an  we  have  already 
shown,  the  value  of  n  has  been  the  prime  factor  in  the  matrix  crack  develop¬ 
ment  processes  observed  in  these  laminates.  Clearly,  it  is  important  to  know 
when  and  how  a  certain  matrix  cracking  pattern  is  reached  at  high  loads,  and 


why  it  is  critical  to  the  final  failure  of  the  load-carrying  plies. 

2.  In  practice,  however,  the  analysts  have  generally  ignored  the  detailed  matrix 
crack  patterns  in  the  off-axis  plies,  but  simply  assume  that  all  loads  will 
eventually  be  taken  by  the  load-carrying  plies.  This  approach  is  typified 
by  the  so-called  progressive  ply  failure  analysis.  The  approach  is  simple 
to  use  and  may  be  adequate  for  preliminary  static  strength  prediction.  The 
reason  may  be  that  matrix  cracks  can  branch  into  the  load-carrying  plies 
only  at  very  high  loads  (except  for  some  extreme  but  impractical  laminates) * 
the  threshold  condition  for  such  event  is  almost  synonymous  to  the  breaking 
of  the  load-carrying  plies  alone.  But,  the  situation  will  be  quite  different 
when  the  applied  load  is  cyclic  in  nature.  Then,  matrix  cracks  in  the  load¬ 
carrying  plies  can  be  developed  at  low  stress  amplitudes  When  this 

happens,  both  the  laminate  residual  strangth  and  fatigue  life  may  not  be 
predicted  based  on  the  properties  of  the  load-carrying  plies  alone.  Then, 
the  detailed  matrix  cracking  patterns  could  become  important*  and  an  accurate 
account  for  their  effects  is  essential.  As  has  recently  been  remarked  by 
Reifsnider,  et.  al.  L5l.  we  quote*  1 

in  order  to  make  an  absolute  (meaning,  quantitative) 
assocaltlon  between  these  (matrix  cracks)  patterns  and 
laminate  residual  strength,  an  accurate  stress  analysis  of 
such  a  pattern  must  be  obtained.  Such  an  analysis  is  Indeed 
complex.  However,  the  nearly  ubiquitous  nature  of  these  types 
of  pattern  in  many  situations  involving  the  long-term  cyclic 
loading  of  composite  laminates  suggests  that  further  study 
of  such  patterns  by  investigators  concerned  with  long-term 
strength  and  life  is  merited." 

This  paper  demonstrated  that  such  an  analysis  is  indeed  complex,  but 


also  sesns  possible 
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SI  Unit 


11 sh  Unit 


144.8  Gpa 

21x10  psi 

11.7  Gpa 

1.7xl06  psl 

0.3 

0.3 

0.54 

0.54 

6.5  Gpa 

0.94x10^  psi 

3.5  Gpa 

0.5xl06  psi 

28.8x10"6/°G 

16x10 “6/°F 

0. 36x10 "V’c 

0.2x10_6/°F 

t,  nominal  ply  thickness 

0.132  mm 

0.0052" 

AT,  temperature  drop 

125  °C 

225  °F 

GIc  (0/0  interface) 

158  j/m2 

0.9  in-lb/in 

GIc  (90/90  interface) 

228  J/m2 

1.3  in-lb/in 

g/t.ttv  gtt/gt=  °*5 

288.75  J/m2 

I.65  ln-lb/i 

See  discussions  In  Ref.  Cl7j 


TABLE  II 


Summary  of  Results  in  Section  2 .1 

Laminate  Type 

1st  Cracking  Event 
(threshold  strain) 

Experimental  value 
(strain  ranges) 

C±25/90]g 

Free  edge  Delam. 

(0.60$ 

Free  edge  Delam. 
(0.58  to  0.61$ 

[02/90]S 

Trans v.  Cracking 
(0.59%) 

Transv.  Cracking 
(  .55  to  0.57%) 

C°2/9o2]s 

Trans v.  Cracking 
(0.33%) 

Transv.  Cracking 
(0.29  to  0.36%) 

C02/90„]3 

Transv.  Cracking 

(0.31%) 

Transv.  Cracking 
(0.24  to  0.28$ 

TABLE  III 


!0~4  IN-LB/IN^r 


Release  Rate  Coeeficient,  C  ,  for  Transverse  Cracking  (AT  =  1  F) 


wlr6oi  £Ni/gi-Ni9oi 


jjLNl/gi-NI  p-01 


Strain  Energy  Release  Rate  Coefficient,  C  ,  for  Transverse  Cracking  (AT  =  1°F) 


MONTE-CARLO  RESULTS  FOR  (0,/90  )c  LAMINATES 


Figure  11  Simulated  Transverse  Crack  Density  vs.  Applied  Tension  Relatioship  for  [o  /90]  . 
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MONTE-CARLO  RESULTS  FOR  (02/902),  LAMINATES 


Simulated  Transverse  Crack  Density  vs.  Applied  Tension  Relationship  for  [02/902]  . 


Simulated  Transverse  Crack  Density  vs.  Applied  Tension  Relationship 


interface 


Lai  Release  Sequence  (below)  for  Delamination  Simulation. 


I06  in-lb/in 


Release  Rate  Coefficient,  C  ,  for  Delsunination  at  Corner  Point 


for  Delami nation  at  Corner  Point. 


Nl/gl-NI  Q 01 


Release  Rate  Coefficient,  C  ,  for  Delamination  at  Corner  Point 


Figure  40  Schematic  View  of  Cross  Cracks  and  Delamination 
at  the  Cross  Cracks. 
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Distribution  on  the  0/90  Interface  Near  the  Cross  Cracks 


I06IN-LB/IN 3 


Strain  Energy  Release  Rate  Coefficient,  C  ,  for  l>el  ami  nation  at  Cross  Cracks 


Strain  Energy  Release  Rate  Coefficient,  C  ,  for  Delamination  at  Cross  Cracks 
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Strain  Energy  Release  Rate  Coefficient,  C-,  for  Delamination  at  Cross  Cracks 


Sequence  of  X  -  Radiographs  Showing  Delami nation  Growth  at  Several  Cross  Cracks 
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